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Abstract 



We present a method for constructing classically perfect anisotropic 
actions for SU(3) gauge theory based on an isotropic Fixed Point 
Action. The action is parametrised using smeared ("fat") links. 
The construction is done explicitly for anisotropy ^ — as /at — 2 
and 4. The corresponding renormalised anisotropies are determined 
using the torelon dispersion relation. The renormalisation of the 
anisotropy is small and the parametrisation describes the true action 
well. Quantities such as the static quark-antiquark potential, the 
critical temperature of the deconfining phase transition and the low- 
lying glueball spectrum are measured on lattices with anisotropy 
^ = 2. The mass of the scalar 0++ glueball is determined to be 
1580(60) MeV, while the tensor 2++ glueball is at 2430(60) MeV. 



^Work supported in part by Schweizerischer Nationalfonds 



1 Introduction 



In Lattice QCD the energy of a physical state is measured by studying the decay 
of correlators of operators which have a non-zero transition matrix element 
between the vacuum and the state under consideration. If the state is heavy, 
these correlators decay very fast in Euclidean time and one has to make sure 
that the temporal lattice spacing is small enough such that the signal of the 
correlator can be accurately traced over a few time slices before it disappears 
in the statistical noise. At the same time, one has to pay attention to choose 
the physical volume of the lattice large enough so that there are no significant 
finite-size effects. Both these requirements together lead to lattices with a large 
number of lattice sites and hence to computationally very expensive simulations. 
The obvious way out of this dilemma is to use a smaller lattice spacing in 
temporal direction compared to the spatial directions, i.e. the use of anisotropic 
lattices. 

Anisotropic lattices have been widely used over the last few years. Studies 
comprising excited states of nucleons heavy-quark bound states ^, 
heavy-light mesons [|[ |^, heavy meson semi-leptonic decays [R, long range 
properties and excited states of the quark-antiquark potentia l [p[ |9|, |ll| as 
well as states composed purely of gluons (glueballs) ^ |14|, [ist or gluons 
and quarks (hybrids) |l^, |9[ |2^, have been performed, mainly us- 

ing the standard anisotropic Wilson discretisation or a mean-link and Symanzik 
improved anisotropic action [ p2[ . In fact, improved actions turned out to be cru- 
cial in anisotropic lattice calculations due to the coarse spatial lattice spacings 
employed. The discretisation of the continuum action, however, can be done in 
many different ways leading to different improvement schemes. A radical ap- 
proach of improving lattice actions, based on Wilson's Renormalisation Group 
ideas, has been suggested by Hasenfratz and Niedermayer namely the use 
of actions that are classically perfect, i.e. there are no lattice artifacts on the 
solutions of the lattice equations of motion. For SU(3) gauge theory the classi- 



cally perfect FP action has been constructed and tested in p4, 25, Eq, E 



and the ansatz has been extended to include FP actions for fermions as well 
In the case of SU(2) gauge theory the FP action has been 
33, pq| and its classical properties have been tested on clas- 



^, ^, 

constructed in 

sical instanton solutions, both in SU(2) and SU(3) |37 . A classically perfect 



gauge action on anisotropic lattices, however, has so far been absent. 

It is thus the goal of this work to fill the gap and to present the construction 
of classically perfect gauge actions on anisotropic lattices. The method has 
been introduced in and it is presented in more detail in |39| . It relies upon 
a recent parametrisation of the isotropic classically perfect Fixed Point (FP) 
action |4^, ^ that includes a rich structure of operators as it is based 
on plaquettes built from simple gauge links as well as from smeared ("fat") 
links. We explicitly construct the ^ — 2 and 4 actions on coarse configurations 
typically occurring in Monte Carlo (MC) simulations. The properties of the 
^ = 2 action are studied by performing measurements of the torelon dispersion 
relation (which serves for determining the renormalisation of the bare (input) 
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anisotropy), of the static quark-antiquark potential, of the deconfining phase 
transition and finally of the spectrum of low- lying glueballs in pure gauge theory. 

It turns out that the construction of anisotropic classically perfect gauge 
actions is feasible. Measuring the renormalised anisotropy using the torelon 
dispersion relation turns out to be stable and unambiguous and shows that the 
renormalisation of the anisotropy is small and under good control. The mea- 
surements of the static quark-antiquark potential indicate that the violations of 
rotational symmetry are small if the (spatial) lattice is not exceptionally coarse. 
This shows that the parametrisation describes accurately the true classically 
perfect action, which is known to have good properties concerning rotational 
symmetry [|^. The study of the glueball spectrum is facilitated a lot due to 
the anisotropic nature of the action, even for (rather small) anisotropy ^ = 2. 
Results, including continuum extrapolations, are obtained for glueball states 
having much larger mass than the highest-lying states that could have been 
resolved with the same amount of computational work using the isotropic FP 
action. However, with the statistics reached so far, it is not yet possible to con- 
clude whether the scaling properties of the glueball states, i.e. the behaviour of 
the measured energies as the lattice spacing is changed, are definitely better for 
the classically perfect action as compared to the mean-link and Symanzik im- 

13| . Especially, this is the case for the lowest-lying 
scalar glueball where the presence of a critical end point of a line of phase tran- 
sitions in the fundamental-adjoint coupling plane, near the fundamental axis, 
causes large distortions of the scaling behaviour (sometimes called the "scalar 
dip"). Our action includes in its rich structure operators transforming according 
to the adjoint representation. If their coupling (which we do not control specif- 
ically during the construction and parametrisation) lies in a certain region, the 
effect of the critical end-point on scalar quantities at certain lattice spacings 
may even be enhanced compared to other (more standard) discretisations with 
purely fundamental operators. Concerning the lattice artifacts observed in the 
glueball simulations, we have to add that other effects, e.g. due to the finite size 
of the lattices, may also be present and that this issue requires further study. 
Being the first application of the FP action technique to anisotropic lattices, 
however, it goes beyond the scope of this work to systematically check all the 
possible sources of errors in the glueball mass determinations. Furthermore, we 
note that the computational overhead of the parametrised anisotropic classi- 
cally perfect action compared to the standard Wilson action as well as to the 
mean-link and Symanzik improved action is considerable. 

The paper is organised as follows. In section ^ we describe our method 
of generating anisotropic perfect gauge actions without having to perform a 
certain number of renormalisation group transformations (RGT) leading from 
very fine to coarse lattices typical for MC simulations. Instead, we make use 
of the parametrisation of the isotropic FP action presented in We 
start with the isotropic parametrised FP action on rather large fluctuations and 
perform a small number of purely spatial blocking steps. If a scale 2 block 
transformation is used, the resulting actions have anisotropics ^ = 2, 4, 8, .... 
The parametrisation ansatz used for these actions is a generalisation of the 



proved anisotropic action 1 12, 
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isotropic parametrisation and is described as well. 

Section ^ contains the main body of the present work, that is, the explicit 
construction and the application of a classically perfect gauge action on ^ = 2 
anisotropic lattices. The results from measurements of physical quantities are 
presented there: the anisotropy of the action is measured using the torelon dis- 
persion relation, the (spatial) scale is determined by the static quark-antiquark 
potential, where for some values of the coupling /?, off-axis separations of the 
quarks are considered as well, in order to estimate violations of rotational sym- 
metry by the parametrised action. Another way of determining the (temporal) 
scale is the study of the deconfining phase transition, measuring the critical 
couplings /?crit for different temporal extensions Nt of the lattice. Finally, the 
glueball spectrum is determined by performing simulations at three different lat- 
tice spacings using the = 2 action and taking the continuum limits whenever 
it is possible. 

The feasibility of iterating the procedure, i.e. adding one more spatial block- 
ing step which leads to a classically perfect action with anisotropy ^ = 4, is 
briefly checked in section ^. It contains the explicit construction and parametri- 
sation of the ^ = 4 action and the determination of the renormalised anisotropy 
via the torelon dispersion relation. 

Section ^ finally contains a summary of the main results, some conclusions 
and an outlook to future work. 

The parameters of the classically perfect actions used throughout this work 
are collected in appendix while appendix |b| contains the detailed results from 
our simulations. 

2 The Construction of Perfect Anisotropic Ac- 
tions 

2.1 Introduction 

The main idea for the construction of (effective) actions on anisotropic lattices 
is to start from an isotropic lattice and to perform one or more purely spatial 
renormalisation group transformations (RGT) . Each scale two RGT step doubles 
the spatial lattice spacing but preserves the spacing in temporal direction. We 
are therefore naturally lead from an action on an isotropic lattice to (effective) 
actions on anisotropic lattices with anisotropics ^ = 2, 4, 8, . . . . Of course, 
choosing the scale s of the RGT different from two, we can in principle reach 
any desired anisotropy ^ = s, s^, s"^, . . . . 

A qualitative picture of the construction of classically perfect actions on 
anisotropic lattices is the following. We start from an isotropic FP action de- 
fined as a point in the infinite dimensional space of all couplings on the critical 
surface where the correlation length is infinite (f3 —^ oo). By performing one (or 
more) spatial RGT we move away from that point always staying on the critical 
surface. Since the spatial blocking kernel is zero on all classical configurations, 
the classical properties of the action are preserved exactly, i.e. it has no lattice 
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artifacts on the solutions of the (anisotropic) lattice equations of motion. (One 
can say that the action is an on-shell tree-level Symanzik improved action to 
all orders in a.) It is, however, important to note that the resulting anisotropic 
classically perfect (AICP) action is, generally speaking, no longer a FP action, 
i.e. the fixed point of a renormalization group transformation. However, it still 
possesses all the good properties of a classically perfect action, in particular it 
preserves scale invariant solutions of the classical equations of motion, i.e. clas- 
sical instanton solutions. Moreover, since the spatial RGT steps are performed 
on the critical surface at /3 ^ oo, the transformation reduces to a saddle point 
problem representing an implicit equation for the anisotropic classically perfect 
action similar to the FP equation in the isotropic case. 

In order to check whether this qualitative picture is correct and that the 
construction works in practice, we have analytically calculated the AICP actions 
for d = 2 scalar field theory (for earlier studies on this subject, see [^) as well as 
for the quadratic approximation to = 4 gauge theory for generic anisotropics 
^. In these two cases it turns out that starting from an anisotropic action for 
very smooth configurations and performing isotropic RGTs or starting from an 
isotropic FP action and performing an anisotropic RGT at the end is equivalent 
and leads to the same classically perfect anisotropic action. It also turns out 
that the quality of the results are as good as what has been obtained with the 
isotropic FP action ^ . In particular we do not see any deviation from the 
classical continuum dispersion relation and almost no violation of the rotational 
symmetry in the pcrturbative potential for r/a > 1. For further details we refer 
toU. ^ ^ 

In practice the main advantage of applying one (or more) anisotropic spatial 
RGT is that we can directly start from a recent parametrisation of the isotropic 
FP gauge action |2^, Q which is valid on rather coarse lattices. We can 
therefore avoid the cumbersome cascade procedure leading from very fine to 
coarser and coarser lattices. 

In the following we first describe the anisotropic spatial blocking in section 



2.2 and then present an extension of the isotropic parametrisation to include 



anisotropy in section 



2.3 



2.2 The Spatial Blocking 

The isotropic FP action is defined by the saddle point equation (at (3 = oo): 

^^P(y) = min {^^P(C/) + TiU, V)} , (1) 

where T{U, V) is the blocking kernel connecting the fine gauge configurations 
U to the coarse configurations V: 

nU,V) = J2 (-^r("s) - ^Re TT[V^inB)Ql{nB;U)]) , (2) 
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where the normalisation term A/'^(nB) is given by 

■^'^'^("^^ ^ vy.Tu^N) {^^' Tr^gt (n^; [/)]} , (3) 

while Q^{tib; U) is the blocked link, 

Qf^ing; U) = W,,(2nB; C/)W^^(2ns + fi;U). (4) 

Here W^{n]U) denotes a smeared (fuzzy) fine link, which is constructed as 
a sum of simple staples as well as of diagonal staples along the planar and 
spatial diagonal directions orthogonal to the link (RGT 3 transformation) |^ , 
respectively. 

To perform the blocking only spatially, doubling the lattice spacing in spatial 
direction while leaving the temporal lattice spacing unchanged, we use a purely 
spatial blocking kernel Tsp{U, V) obtained by setting in the isotropic definition 

Qi{nB;U) = Wii2nB;U) (5) 

thus doing only a smearing and no blocking in temporal direction. (In the case 
of spatial blocking the expression "2^3" stands for {2n^ , 2n^ , 2n^ , n%) .) In 
addition, we have the freedom to choose different values of k for the spatial (ks) 
and temporal (nt) blocking. 



2.3 The Parametrisation 

The classically perfect actions, described in principle by an infinite number 
of couplings, have to be parametrised in order to be suitable for numerical 
simulations. We have shown previously j28[ ^ that a parametrisation of 
the isotropic FP gauge action including APE-like smearing behaves much better 
and is much more flexible compared to common parametrisations using traces 
of closed loops with comparable computational cost. We thus decide to use an 
extension of this parametrisation to describe the anisotropic classically perfect 
gauge action. 

For convenience, let us quickly review the parametrisation of the isotropic 
action. The parametrisation uses mixed polynomials of traces of simple loops 
(plaquettes) built from single gauge links U^{n) as well as from asymmetrically 
(APE-like) smeared links w'l^\n). We introduce the notation sj^\n) for the 
sum of two staples of gauge links connecting two lattice sites in direction ^ lying 
in the fiv-jAane: 

^ U,{n)Uf,{n + v)Ul{n + ii) + Ul{n~ i>)Uf,{n - v)U,{n - v + (i). 

(6) 

To build a plaquette in a plane from smeared links it is convenient to 
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introduce asymmetrically smeared links. First definqj 

Q'^^'-i[t. Sl''+vSlr'y{^ + lv)u,. (7) 

Out of these sums of matrices connecting two neighboring points rt, n + fl, we 
build the asymmetrically smeared links 



i=l 



The parameters rj, Ci used for the smearing may depend on local fluctuations 
measured by x^(n) defined as 

x^{n) = Re Tr(Q^(n)C/t(n)), (9) 

with the symmetrically smeared link 

Q:<N = ^E^m'H^)-c^a.(")- (10) 

This parameter is negative, —4.5 < Xfj_{n) < 0, and it vanishes for trivial gauge 
configurations. 

The smearing parameters are chosen to be polynomials of with free coef- 
ficients (determined later by a fit to the FP action): 

77 = 77(0) +77(1)2; + 77(^)2:2 +... , (11) 

(12) 



Of course, the asymmetrically smeared links VF/i"^^ built out of a large number 
of paths connecting the neighbouring lattice sites are no longer elements of the 
SU(3) gauge group. They might be projected back to SU(3), however this task 
increases the computational cost in actual numerical simulations. Moreover, our 
studies have shown that projection reduces the degrees of freedom in defining 
the action and we are thus using the smeared links W^i'^'' as they are. 

We now build a smeared plaquette variable, 

w^, = ReTT{l-WPl), (13) 

as well as the ordinary one, 

w^, = ReTr(l-f/Pi), (14) 

where 

WPlin) = VK^'^'(7i)W^^^'(77 + A)W^i''^^(" + i>)W^i^'^(n) (15) 



The argument n is suppressed in the following 
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and 



U^lin) - U^{n)UUn + A)f/^(" + i^WHn), (16) 

respectively. Finally, the action is built from these plaquette variables using a 
mixed polynomial ansatz of the form 

^[^] = ]^ E E E^'fc'"''''H'^M^(")'' (17) 

n ^i<u k.l 

where the coefficients pki are again free parameters defined by a fit to the FP 
action. 

In order to adapt this parametrisation of the isotropic FP action to anisotropic 
actions, we use three (rather straightforward) extensions. Firstly, the coeffi- 
cients pki in eq. ( p^ ) are chosen differently depending on the orientation of the 
plaquette fxv, i.e., for fjLv S {12,13,23} (spatial plaquettes) and for 
jiv G {14,24,34} (temporal plaquettes). Secondly, the parameter 77 entering in 
eq. (0) describing the asymmetry between differently oriented staples contribut- 
ing to a smeared link wj^ ^ is generalised to distinguish between smeared spatial 
and temporal links: 



l(4) 



\{ E ' + ^1^1'^ + ^35f ^) - ^(1 + ^1 + m)u., (18) 

i(E #'+'?4^fV(l + ^'?4)t/., (19) 

^4^ - \(Y.st^ + ms'i^)^{i + \m)u,, (20) 

where i, j, k — 1,2,3 and — \^ . . . ,4. The anisotropic parameters rji, . . . ,7^4 

may be again polynomials in the local fluctuation parameter a;^. These situa- 
tions are depicted in Figure |l]. 

Finally, the construction of the smeared links wjf^ from the matrices qI^"* , 
described by the parameters Ci in eq. (H), is generalised such that these pa- 
rameters are chosen differently for the construction of temporal links (always 
contributing to temporal plaquettes), spatial links contributing to spatial pla- 
quettes and spatial links contributing to temporal plaquettes, qi, Ci2 and 0^3, 
respectively, see Figure ^ 

It is not a priori clear whether all the three extensions are necessary at the 
same time. Results from fits to different extended parametrisations (Table 
however, show that it is indispensable to have the full parameter set. The values 
of of the fit to the true classically perfect actions as well as t he measured 



renormalised anisotropics of the parametrised action (see Section 3^) indicate 
that the parametrisation is not flexible enough neither for 773 = 774 = (leading 
to a positive definite transfer matrix connecting only neighbouring time slices) 
nor if the parameters c^j, j = 1,2,3 are set to be equal. 
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m 



V3 



Figure 1: Asymmetry in the construction of the smeared hnk matrices Q^T''- 
From left to right: spatial staple contributing to a smeared spatial link in a 
spatial plaquette (?7i), temporal staple contributing to a smeared temporal link 

(772) , temporal staple contributing to a smeared spatial link in a spatial plaquette 

(773) , temporal staple contributing to a smeared spatial link (774). 



Cil 



Ci2 



J 

Ci3 



Figure 2: Asymmetry in the construction of the smeared links W^^' from the 
matrices Q^T^- From left to right: spatial Qp-* contributing to a smeared spatial 
plaquette {cn), spatial Qf^ contributing to a smeared temporal plaquette (0^2), 
temporal Q'f^ (^3). 



Our anisotropic action shall have the correct normalisation and anisotropy 
in the continuum limit. This can be accomplished by demanding the following 



two normalisation conditions to be exactly fulfilled |39 



Po\+Pfo + 2PoiCl3+PoiCl3?72 + 2p^lCll??3+PoiCl2'74 = (21) 

P"i+/io + 2poiCii + 2p^\ci2 + 2p^^iCii77i = (22) 

where all parameters 77, c denote the constant (O"^ order) terms in the polyno- 
mials in x^. 

Let us note that the parametrisation presented here leads to an overhead in 
the updates in MC simulations of 66 compared to the Wilson action and of 12 
compared to the tadpole and tree-level improved action. 



3 The ^ =2 Perfect Action 
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#c. 


n 


inax(fc + l)sp 


niax(A; + Z)tm 


9 




4 


1 


4 


4 


4 


0.0250 


1.63(2) 


2 


3 


3 


4 


4 


0.0238 




4 


3 


3 


4 


4 


0.0144 


1.912(9) 



Table 1: Comparison of the accuracy Xd of the fit on derivatives on 20 con- 
figurations at /3 = 3.5 and the measured renormahsed anisotropy (using 
the torelon dispersion relation) at /3 = 3.3 for different choices of the set of 
non-linear parameters. 



3.1 Construction 

We construct the ^ = 2 classically perfect anisotropic action starting from 
the parametrised isotropic action HI] performing one spatial blocking step 
as described in Section ^. Studies of the quadratic approximation suggest 
Kt = £.^Ks — 4ks in the RG transformation for keeping the block transfor- 
mation as close to the isotropic case as possible. Furthermore, spatial locality 
should stay at its optimum for the "isotropic" value of Ks = 8.8 ||2^, thus we 
choose Kt = 35.2. Indeed, it turns out that varying Kt away from Kt — 35.2, 
keeping Ks = 8.8 fixed, makes it much more difficult to parametrise the resulting 
perfect action such that the renormalised anisotropy stays close to the input 
value ^ = 2. Coarse configurations V are generated using the ad- hoc anisotropic 



action (see Appendix A. 2), the FP eq. (M including the purely spatial blocking 



kernel Tsp{U,V) is used as a recursion relation with the intermediate isotropic 



action (see Appendix A.l) on the r.h.s. yielding the classically perfect = 2 ac- 
tion on the l.h.s. Examining the blocking, we observe that coarse configurations 
obtained in a MC run with the ad-hoc anisotropic action are mapped to fine 
configurations that are close to isotropic (concerning spatial and temporal pla- 
quette values and expectation values of the Landau gauge-fixed link variables) 
for input anisotropy ^ad-hoc ~ 3.2. We thus generate 20 configurations each at 
/3ad-hoc = 2.5, 3.0, 3.5, 4.0 using the ad-hoc anisotropic action^with ^ad-hoc=3.2. 
These configuratio ns ar e spatially blocked using the intermediate isotropic ac- 



tion (see Appendix A4 ) which describes the minimised isotropic configurations 
on the r.h.s. of eq. (|l|) reasonably well. 

To construct the action, we perform several non-linear fits for the derivatives 
of 20 configurations at /3ad-hoc = 3.5, using different sets of parameters. The 
parameters Ci are chosen to be non-zero for j = 1, 2, 3 as adding additional free 
C4 parameters does not improve significantly. 

^Note, that the values of /3ad-hoci used in this section, may not be directly compared to 
the values of /3perf corresponding to the parametrised classically perfect action, used in the 
following sections; but rather the fluctuations (plaquette values, expectation values of Landau 
gauge-fixed links etc.) of the gauge configurations produced in MC runs using the perfect 
action should lie approximately in the same range as the fluctuations of the initial coarse 
configurations produced with the ad-hoc action. 
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Having fixed the non-linear parameters rj and we include the derivatives 
of 20 configurations at /3ad-hoc = 4.0 and the action values of all the configura- 
tions at /3ad-hoc = 3.5, 4.0 in the linear fit of the parameters pki, keeping the 
non-linear parameters fixed. The resulting values of and the linear behaviour 
of the actions suggest a linear set with max(fc-|-Z)sp = 4, max(fc-|-Z)tm = 3 where 
the action values are included with weights Wact = 0.018 relative to the deriva- 
tives. Using this value, the good parametrisation of the derivatives obtained in 
the full nonlinear fit is preserved, at the same time the mean error of the action 
value due to the parametrisation is as small as 0.35%. It is checked that the 
number of data points in the fit is large enough, making sure that the values 
of as well as the errors of the action values do not increase significantly on 
configurations independent of those used in the fit. Furthermore, these quan- 
tities are determined on configurations down to /3ad-hoc = 3.0 and it turns out 
that the action works down to this value of /3ad-hoc. Finally, it is checked that 
the linear parameters pki do not lead to artificial gauge configurations in MC 
simulations caused by (local) negative contributions to the total action. The 



parameters of the resulting action are listed in Appendix A. 3 



3.2 The Renormalised Anisotropy 

The renormalised anisotropy of the action presented above is measured using 
the torelon dispersion relation as described in The torelon is a closed 

gluon flux tube encircling the periodic spatial boundary of size L of a torus. 
It is created and annihilated by closed gauge operators winding around this 
boundary. String models which are a good approximation for long flux tubes 
predict its energy to behave as 

E{L)^aL + ^, (23) 

where a is the string tension and the second term describes the string fluctuation 
for a bosonic string |6[ |4^, The string models assume no restriction 
on the transverse modes of the string, thus they apply to volumes with large 
transverse extension S. Due to this reason and because we prefer not to have 
too large transverse momenta, we work on S'^ x L x T lattices, where S and T 
are large, whereas the length of the z-dircction L is chosen to be rather short 
such that the torelon energy (see eq. ([23|)) is not too heavy. Since the colour 
flux tube has finite width in physical units, it is useful to employ (iteratively) 
APE smeared links to improve the overlap of the operators with the ground 
state. One step of APE smearing acts on a spatial gauge link Uj{x) as 

SiUj{x) = Vsv(^s^[u,{x)+\sY,iUk{x)U,{x + k)ul{x + j) (24) 

+Ul{x - k)U,{x - k)Ukix -k + j))}, 

the original link variable is replaced by itself plus a sum of the four neighbouring 
spatial staples and then projected back to SU(3). The smeared and projected 
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links SiUj{x) have the same symmetry properties under gauge transformations, 
charge conjugation and reflections and permutations of the coordinate axes as 
the original gauge links. The smearing is used iteratively allowing to measure 
operators on different smearing levels, SnU . 
We measure the operators 

L-l 

Tn{x,y,t) ^TvW_SnU:,{x,y,z,t), (25) 

z=0 

for all X, y, t, where SnUz is the n times iteratively smeared link in the longitu- 
dinal z direction. By discrete Fourier transform, we project out the state with 
momentum p = {px,Py) = {nx,ny)(2n/Sas): 

T„(p, t)^Y. ^"(^' 2/' t)e'^^^"+''""^ (26) 

and build the correlators 

C„„,(p,t) - ^^(r„(p,T)T„,(-p,T + i)). (27) 

From these correlators of operators measured on the different smearing levels 
we obtain the energy values atE{p) using variational methods (see Appendix B 
of [^sl ). The continuum dispersion relation 

= p^ + (28) 

on the lattice becomes (in temporal units) 

(29) 



where fix 5 fly Q-rG the components of the (transversal) lattice momentum. On an 
anisotropic lattice, this equation allows for the extraction of the renormalised 
anisotropy — as /at (measuring the "renormalisation of the speed of light") 
as well as the torelon mass mof, which in turn may be used to get an estimate 
of the scale using eq. (|2^) and known values of the string tension a. The main 
advantage of this approach is obtaining the renormalised anisotropy ^j? as well 
as an estimate for the scales and at by performing computationally rather 
inexpensive measurements. 

The parameters of the simulations carried out are collected in Table |^. 

The Polyakov line around the short spatial direction is measured using 
operators at smearing levels 5„C/ with n = 3,6, 9, 12, 15 for (3 > 3.3 and 
n = 2,4,6,8,10 for fi < 3.15, respectively, and smearing parameter As — 0.1. 
The measured energies Eljp') are given, together with the number of operators 
used in the variational method, the fit ranges and the values of per degree 
of freedom, x^/^df, in Tables p3 and in Appendix M. 
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p 


S'^ X LxT 


# sweeps / measurements 


3.00 


8^ X 4 X 20 


72000 / 14400 


3.15° 


8^ X 5 X 20 


54000 / 10800 




52 X 8 X 20 


47600 / 9520 


3.30 


12^ X 6 X 30 


35600 / 7120 


3.50° 


12^ X 8 X 24 


39400 / 7880 


3.50'' 


142 X 6 X 30 


36800 / 7360 



Table 2: Run parameters for the torelon measurements using the ^ = 2 perfect 
action. The lattice extensions in torelon direction L, the extension in the two 
transversal spatial directions S as well as the temporal extension L are given. 



Figure ^ displays an example of a dispersion relation at (3 = 3.3 including 
all values of p at which the energies could be determined on the given lattice. 
Note, that the energies of the p — torelons (the torelon masses) may be hard 
to determine because the effective masses do not reach a plateau within the 
temporal extent of the lattice, probably due to the small overlap of the operators 
used with these states; however, using the non-zero momentum energies, the 
masses may still be accurately determined. 

To determine the renormalised anisotropy as well as the torelon mass 
in units of the temporal lattice spacing mat, we perform fits to the dispersion 
relation, taking into account the correlations between different operators. The 
range of considered is chosen depending on x^/-^df of the fit and the precision 
of the dispersion relation data in the respective range. The results are given in 



Table |. 


f3 


fit range 




rriTat 




3.00 


0..5 


1.903(81) 


1.324(98) 


0.39 


3.15° 


0..4 


1.966(39) 


0.700(35) 


0.83 


3.15'' 


0..5 


2.022(104) 


1.262(84) 


0.62 


3.30 


1..9 


1.912(9) 


0.311(5) 


0.94 


3.50° 


1..10 


1.836(9) 


0.149(10) 


0.94 


3.50'' 


1..8 


1.826(16) 


0.208(16) 


1.61 


Table 3 


: Results of the torelon simulations 


using the ^ = 


2 perfect action. The 


fit rang 


e in p^ is given 


in units of (2TT/Sf. 
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Using the finite size relation for the torelon mass corresponding to eq. , 

MriLa,) = (a + -J^)Las , (30) 
(Las ) 

we calculate the string tension a and thus the spatial scale as ■ The string model 
predicts D = — ^ for long strings thence we use this string picture value and 
stay aware that the estimate gets worse for short strings. To obtain the hadronic 
scale ro and the lattice spacings we employ ro-\/a'=l-193(10) from |^ and use 
the definition tq = 0.50 fm. The results are collected in Table ^ 

The two simulations denoted by 3.15'' and 3.50'' are used to check the stabil- 
ity of the method. The first one is carried out on a lattice with small transversal 
spatial size (the torelon string measured is even longer), the latter is measuring 
the energies of a rather short string. The scale is measured accurately using 



the static quark-antiquark potential (see Section 3.3.1) so that we can compare 
our torelon estimates to the much more reliable values in Table ^. It turns 
out that the deviation of the torelon estimates for the scale ro /as from the po- 
tential values does not exceed 6% if the length of the string is not too short 
(say Las ^ 1 fm). For short strings, the string picture does not apply and 
we may therefore not expect good results in this case. But even for very long 
strings, it is not possible to obtain very accurate information about the scale 
as long strings are very heavy and thus difficult to measure (comparable to the 
long range region of the static quark-antiquark potential), especially when the 
anisotropy ^ is not very large. 

There seems to be no problem if the transverse volume is rather small as in 
the (3 — 3.15'' simulation, except of course the large momenta occurring which 
make the determination of the energies more difficult. As well, the estimates 
of the renormalised anisotropy do not show significant deviations neither for 
small transverse volume nor for short strings. 



p 


y/aas 


ro/as 


as [fm] 


at [fm] 


3.00 


0.834(20) 


1.43(5) 


0.350(12) 


0.184(14) 


3.15" 


0.563(8) 


2.12(5) 


0.236(6) 


0.120(5) 


3.15'' 


0.579(18) 


2.06(8) 


0.243(9) 


0.120(11) 


3.30 


0.3578(21) 


3.33(5) 


0.150(2) 


0.079(1) 


3.50'^ 


0.225(5) 


5.31(16) 


0.094(3) 


0.051(2) 


3.50'' 


0.304(7) 


3.92(12) 


0.128(4) 


0.070(3) 



Table 4: Estimates of the scale determined from torelon results (see Table 
see text. 



Finally, we may conclude that determining the renormalised anisotropy using 
the torelon dispersion relation is a stable and apparently sensible procedure. 
There are no manifest problems and it works on fine lattices as well as on 
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13 Lag [fm] 


7'o/as (Torelon) 


ro/os (Potential) 


rel. error 


3.00 1.48 


1.43(5) 


1.353(16) 


+5.7% 


3.15° 1.23 


2.12(5) 


2.038(2) 


+4.0% 


3.15'' 1.96 


2.06(8) 


2.038(2) 


+ 1.1% 


3.30 0.95 


3.33(5) 


3.154(8) 


+5.6% 


3.50" 0.82 


5.31(16) 


4.906(21) 


+8.2% 


3.50'' 0.61 


3.92(12) 


4.906(21) 


-20.0% 



Table 5: The estimates of the scale obtained from torelon measurements (see. 
Table ^ compared to the scale measured using the static quark-antiquark po- 
tential (see Table ||). The length of the torelon string Lag in physical units is 
given as well. 




J I I I I I I I I I I I I I \ I \ I L 

2 4 6 8 10 12 14 16 18 
(a,p)" in units of (27i/12)^ 



Figure 3: Torelon dispersion relation for (3 = 3.3. The straight line is the 
correlated fit to E'^{p) = + in the range — 1..9. 
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coarse ones. Checks employing different lattice extensions indicate that the 
determination of the renormalised anisotropy is quite stable within a large range 
of different choices. The estimate of the lattice scale, however, is not very 
accurate, mainly due to the use of the string picture relation, eq. (^0|); the 
(systematic) error of the torelon scale is about 5% for reasonable torelon lengths 
L > 1 fm with a tendency of underestimating the lattice spacing Qs ■ 

3.3 The static quark-antiquark potential 
3.3.1 The Scale 

To set the scale of the action at a given /3-value we measure the static quark- 
antiquark potential, which is reasonably well described by the phenomenological 
ansatz [|o|, 

rv 

V{R) ^Vn + -+aR. (31) 
K 

We determine the lattice spacing using the following definition of Vc |^ : 

r^V'ir)\r^r^=c, (32) 

where the well-known Sommer scale tq ~ 0.50 fm corresponds to c = 1.65. For 
the anisotropic action, we employ the conventional definition of tq, however, 
V{r) is now measured in units of the temporal lattice spacing and thus we have 

-a + C<T( — )2 =c, (33) 

as 

where a denotes the fitted dimensionless value a = aa^at of the string tension. 
The renormalised anisotropy ^ thus has to be known before the spatial scale 
To may be determined. For very fine or coarse lattices it may be advantageous 
to go to a slightly different separation corresponding to a different value of the 
constant c on the r.h.s. For this purpose, we have collected values for c and 
Tc from high-precision measurements of the static potential performed with the 
Wilson action Q listed in Table ^. 



rc/ro 


c 


0.662(1) 


0.89 


1.00 


1.65 


1.65(1) 


4.00 


2.04(2) 


6.00 



Table 6: Parameter values for the determination of the hadronic scale Tc through 
r^V'ir)\r=r^=c. 

The static quark-antiquark potential is measured on lattices of different 
scales fls = 0.10 . . .0.37 fm, employing Wilson loops built out of APE smeared 
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(see Section |3.2| ) spatial links and simple temporal gauge links. The measure- 
ments at (3 = 3.00, 3.30 include the determination of the potential between 
quarks that are ofF-axis separated along the lattice vectors (1,0,0), (1,1,0), 
(1,1,1), (2,1,0), (2,1,1), (2,2,1) (and lattice rotations thereof) in order to es- 
timate violations of rotational symmetry. The rest of the measurements include 
only on-axis separations due to the large computational cost (concerning speed 
and memory) of the off- axis measurement. The smearing levels used are as fol- 
lows (A/ = 0.1 on all the lattices): on the finest lattice aX (3 — 3.50, we use 
SnU,n = 3,6,9,12,15; at f3 — 3.15 we use SnU,n — 2,4,6,8,10; for the mea- 
surements including off-axis separations, we employ only three different smear- 
ing levels each, in order to save memory and time: SnU, n — 5, 10, 15 at /3 = 3.30 
and SnU, n = 2, 4, 6 on the coarsest lattice at /3 — 3.00. The parameters of the 
simulations are given in Table |. 





off-axis Sep. 


53 X T 


# sweeps / measurements 


3.00 


yes 


8^ X 16 


42000 / 2800 


3.15 


no 


10^ X 20 


39800 / 3980 


3.30 


yes 


10^ X 20 


27000 / 1800 


3.50 


no 


123 X 24 


36400 / 3640 



Table 7: Run parameters for the measurements of the static quark-antiquark 
potential using the ^ = 2 perfect action. 



The values of the off-axis potential atVir") at f3 ~ 3.30 are collected in Table 
p6| in Appendix ^ and displayed in Figure ^, the values of the ofF-axis potential 
on the coarse lattice at (3 = 3.00 are collected in Table p7|. The potential values 
atV{r) of the on-axis simulations are collected in Table |28|. 

The parameters a and a (the string tension) in the phenomenological po- 
tential ansatz, eq. (]3l]), are determined using global fits, the hadronic scale tq is 
determined performing local fits. The fit ranges and results are given in Table 
^. In the global fit to the off-axis potential at P — 3.30 we exclude the badly 
measured separations 2(2, 1, 0) and 2(2, 1, 1). 

In principle, it is possible to determine the renormalised anisotropy using the 
first excited state of the potential together with the ground state. The ground 
state values are again fitted to the ansatz V{r) — Vq + a/r + ar, whereas the 
values of the excited state are fitted to the (ad-hoc) ansatz V*{r) = A/r + B + 
Cr + Dr^. As can be seen from Figure ||, however, for the off-axis potential 
at /3 = 3.3, the energy values of the first excited state have large errors. The 
result for the separation at{V*{ro) — ^(''o)) = 0.555(47) thus shows a large 
error. Comparing this value to ro(y*(ro) — V{ro)) w 3.25(5) of |Q one obtains 
^ = 5.86(59) and thus = Os/at=l. 86(19) for (3 = 3.3 which is in agreement 
with the value determined using the torelon dispersion relation, ^r=1. 912(9). A 
similar determination for the on-axis potential at /? = 3.50 is even more difficult 
due to the smaller number of separations measured. The result is 1.66 (32) 
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0.8 



0.6 - 



L 0.4- 



0.2 




Figure 4: Ground state of the static quark-antiquark potential at /3 = 3.3. 
Including a fit to the phenomenological potential ansatz in the range r — y/2..5. 



p gl. 


fit 


a 




loc. 


fit 


ro/as 




3.00 1 - 


2V6 


-0.201(18) 


0.3595(87) 


V6- 


- 3 


1.353(16) 


1.119(51) 


3.15 1- 


- 5 


-0.1503(12) 


0.1654(6) 


1 - 


3 


2.038(2) 


1.162(15) 


3.30 V2 


- 5 


-0.1539(7) 


0.0683(3) 


2^/2- 


2V3 


3.154(9) 


1.140(8) 


3.50 1- 


- 6 


-0.1478(4) 


0.0320(3) 


2- 


4 


4.906(21) 


1.189(14) 



Table 8: Results of the measurements of the static quark-antiquark potential 
using the ^ = 2 perfect action. The parameters a and a of the potential ansatz 
eq. (|3^) are given together with the global fit range chosen to determine them. 
The hadronic scale rg is given in spatial lattice units together with the local fit 
range used to determine that quantity. Note, that we have chose n c = 6.00 and 
c = 0.89 for /3 = 3.0 and (3 = 3.5, respectively (see Section 3.3.1 ). Additionally, 
the dimensionless quantity r^yfa is given. 
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which does again agree with the torelon result, <^r, = 1.836(9). Due to the large 
errors (possibly caused by the operators used, optimised for the ground state) 
this method of determining the renormalised anisotropy is not suitable in 
the context of this work. However, for the determination of the scale as well as 
of the anisotropy on lattices with finer temporal lattice spacing (above all with 
higher anisotropics) this way might be feasible at least if off-axis separations are 
included in the measurement. 




Figure 5: Static quark antiquark potential at (3 ~ 3.3 including the first excited 
state V*{r) at /3 = 3.3. The location of tq where the determination of the 
anisotropy is done, using the gap, is marked by the dotted line. 



3.3.2 Rotational invariance and scaling 

The measurements at f3 — 3.0, 3.3 including a large number of off-axis separa- 
tions provide information about the deviations from rotational invariance. This 
issue has been addressed in for the RGT used in this work and it turned 
out that violations of rotational symmetry caused by the blocking are small. 
That this is still true for the parametrisation employed here is shown by the 
off-axis potential at (3 = 3.30, see Figure ^. The deviations on the coarse lattice 
at (3 = 3.00 are significantly larger (see Figure ^ which is no surprise as this 
coupling corresponds to a spatial lattice spacing w 0.37 fm. We present some 
dimensionless ratios containing information about the deviation from rotational 
invariance in Tab. ^; due to the resulting small values the errors are rather 
large. One notices that the parametrised perfect action does a good job, for 
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f3 — 3.30 the deviation between the exactly degenerate separations (2,2,1) and 
(3,0,0) is zero (within statistical errors). 



f3 


V(l,l,l)-Vfit(\/3) 

y(2,o,o)-y(i,o,o) 


y(2,2,i)-y(3,o,o) 
y(3,o,o)-y(i,o,o) 


3.00 


0.086(4) 


0.057(6) 


3.30 


0.029(7) 


0.005(6) 



Table 9: The normalised deviations from rotational invariance for the AICP 
action. The ratios displayed are zero in the case of no rotational symmetry 
violation. 



The static gg-potential is also an effective test of scaling. Expressing the 
potential measurements performed at different couplings /3 in the RG invari- 
ant, dimensionless ratios r/rp and rgV and subtracting the unphysical constant 
roV{ro) should lead to potentials lying exactly on top of each other. Deviations 
indicate either scaling violations or ambiguities in the determination of tq. Fig- 
ure ^ includes all potential measurements and shows that the different curves 
can hardly be distinguished from each other, except some energies determined 
on the coarsest lattice at /3 = 3.00. They deviate notably from the curve which 
is a fit to the phenomenological potential ansatz, eq. (^ij), including all the mea- 
surements that have been included into the global fits of the single f3 values, 
see Table ||. The results for the two physical parameters from this global fit are 
a = -0.27799(5) and crrg==l. 3690(3). 



3.4 The Critical Temperature 

The critical temperature of the deconfining phase transition contains informa- 
tion about the temporal scale of the lattice at a given coupling /3 as Tc = 
l/{Ntat). Comparing this information to quantities obtained from the measure- 
ments of torelons or the static quark-antiquark potential offers many interesting 
scaling (and other) tests. However, since measurements of such quantities at 
the determined values of Pent are yet absent and because interpolations in /3 are 
very difficult due to the renormalisation of the anisotropy ^, these tests are a 
future project. 

The critical temperature Tc may be defined as the location of the peak in 
the susceptibility of the Polyakov loop P (Wilson line in temporal direction), 

X^N!{{\Pf)^{\P\f). (34) 

We determine the values of Pc{Nt) on rather coarse lattices with temporal 
extensions A^t = 3, . . . ,7. Simulations at different values of P near the estimated 
critical couplings Pc are performed and the peaks are determined by employing 
the Ferrenberg-Swendsen multi-histogram reweighting js^, |5^ (see also ^ 
^ for additional details). The run parameters of the simulations are collected 
in Table ITol. 
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Figure 6: The static quark- antiquark potential measurements expressed in RG 
invariant, dimensionless units r/ro, r^V. The unphysical constant roV{ro) has 
been subtracted such that aU the curves exactly coincide at r = tq. The dashed 
line is a global fit to the phcnomcnological potential ansatz, cq. (pl|). 
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In this first study, we decide not to examine the finite-size scaling of /Scrit 
but to choose Lg/Lt = ^ ■ Ns/Nt « 3.5 ~ 4. Obviously, the anisotropic nature 
of the action makes the computational effort of the simulations for exploring 
the deconfining phase transition smaller. The effects of the finite volume are 
(rather conservatively) estimated by using the finite-size scaling relation for a 
first order phase transition. 



The parameter \h\ is set to the largest value \h\ = 0.25 that has been observed 
with the isotropic action, and the value of ^r, used to estimate the spatial lattice 
volume is set to the lowest value appearing in the whole range of f3 values 
considered. The extrapolated critical couplings (3c are given together with the 
statistical error and the estimated finite- volume error in the last column of Table 



3.5.1 Introduction 

The particles mediating the strong interaction of QCD, the gluons, carry colour 
charge and thus interact with each other, unlike, e.g., their counterpart in elec- 
tromagnetism, the photons, which have zero electric charge. The spectrum of 
QCD may thus contain bound states of (mainly) gluons, called glueballs. These 
states are described by the quantum numbers J denoting the (integer) spin, P 
denoting the eigenvalue ±1 of the state under parity and C denoting the eigen- 
value ±1 under charge conjugation. Thus the eigenstates of the Hamiltonian 
corresponding to glueball states are labeled | J^*^). 

Currently, there is an ongoing debate whether light glueballs (above all the 
scalar 0++ which is the lightest state in pure lattice gauge theory with a mass 
of about 1.6 GeV) have been observed experimentally at about the mass that is 
predicted by quenched simulations on the lattice, whether the lightest glueball is 
much lighter (below 1 GeV) and very broad |5^, or whether glueballs have not 
been observed at all in experiments ||59|| . There are mainly two reasons for this 
uncertainty. On one hand, the experimental data seem not yet to be accurate 
and complete enough, despite large efforts in the last years, driven by the lattice 
results; on the other hand, lattice simulations with high statistics, measuring 
glueball states, have been performed only in the quenched approximation, where 
the quarks are infinitely heavy and thus static. Decreasing the sea (dynamical) 
quark mass (finally down to the physical value) will allow to track the glueball 
states as sea quark effects are increased. It may turn out, that indeed the 
glueball mass is lighter than the one measured in pure gauge theory (for partially 
quenched results possibly indicating this see ^ ^)- It may even happen 
that by "switching on" the sea quarks the scalar glueball starts to decay (almost) 
instantaneously to qq states, i.e. it ceases to exist physically. However, the 
(partially) unqucnchcd results are rather indecisive yet. 




(35) 
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Glueballs 
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lattice size 





sweeps 


Tn 
P 






0c 


3 X 6^ 


2.80 


45000 




70 


1.45(5) 


2.863(1)(5) 




2.85 


45000 


410 


259 


11.3(20) 






2.87 


45000 




312 


21.6(28) 






2.90 


50000 




100 


2.34(10) 




4 X 8^ 


3.00 


45000 




81 


1.63(5) 


3.032(1)(5) 




3.02 


45000 


870 


102 


3.97(26) 






3.03 


45000 


830 


372 


43.2(22) 






3.04 


45000 




188 


10.9(18) 






3.05 


55600 




144 


3.94(21) 




5 X 93 


3.07 


44000 




117 


0.84(4) 


3.118(1)(6) 




3.10 


45000 


560 


154 


4.38(63) 






3.11 


20500 


640 


198 


15.0(27) 






3.12 


44500 


3700 


278 


33.2(19) 






3.15 


20000 




110 


3.26(18) 




6 X 11^ 


3.17 


45000 


1440 


121 


7.4(11) 


3.181(1)(6) 




3.18 


45000 


1750 


201 


30.0(15) 






3.185 


45000 




232 


23.2(24) 






3.19 


45000 




123 


10.4(14) 






3.20 


45000 




100 


5.62(36) 




7 X 13^ 


3.22 


32400 


700 


87 


6.7(5) 


3.236(1)(6) 




3.23 


44500 


750 


129 


18.4(11) 






3.24 


29000 




149 


20.7(16) 






3.25 


22400 




83 


9.2(10) 





Table 10: Run-time parameters of ^ = 2 Tc-simulations. The number of sweeps, 

the persistence time Tp (if applicable), the integrated autocorrelation time Tint 
and the Polyakov loop susceptibility Xh are given for all values of used in the 
final reweighting procedure, as well as the resulting 0c. The systematic error 
in the second brackets is a (rather conservative) estimate of the finite volume 
effects and has to be considered if the given value should be an estimate of the 
infinite-volume critical coupling. 
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pp 


Do 


Di 


D2 


^3 


Di 


D5 


De 




1 











1 





1 


A2 











1 








1 


E 








1 





1 


1 


1 


Ti 





1 





1 


1 


2 


1 


T2 








1 


1 


1 


1 


2 


Table 11: 


: The composition 


of the subduced r 


epresentation Dj 


iOh 


in terms of 



the irreducible representations of the cubic group Oh- 
3.5.2 Glueballs on the Lattice 

Glueballs in the continuum are rotationally invariant and have a certain (integer) 
spin J. On the lattice, the rotational 0(3) symmetry is broken, only its discrete 
cubic subgroup Oh survives the discretisation. Therefore, the eigenstates of the 
transfer matrix are classified according to the five irreducible representations 
of Oh'- Ai, A2, E, Ti, T2 with dimensions 1, 1, 2, 3, 3 respectively. Their 
transformation properties may be described by polynomials in the components 
X, y, z of an 0(3) vector as follows: Ai ~ {1}, A2 ^ {xyz}, E ^ {x^ — , — 
z^}, T\ ~ z}, T2 ^ {xy, xz,yz}. Generally, an 0(3) representation with 

spin J splits into several representations of the cubic group. Since Oh is a 
subgroup of 0(3), any representation Dj with spin J in the continuum induces 
a so-called subduced representation Dj J, Oh on the lattice. This subduced 
representation no longer has to be irreducible but is a direct sum of irreducible 
representations of Oh- 

Dj iOh^T^(Br^(B--- - (36) 

Table |ll] Hsts the subduced representations of Dj for J — 1, . . . ,6. The spin 
J = 2 state for example splits up into the 2-dimensional representation E and 
the 3-dimensional representation T2. Approaching the continuum, rotational 
symmetry is expected to be restored and thus, as a consequence the mass split- 
ting of these two states will disappear and the two representations form together 
the 5 states of a spin J — 2 object. 

Pure glue physical states on the lattice are created and annihilated by apply- 
ing gauge invariant operators to the pure gauge vacuum. In our simulations, we 
use space-like Wilson loops in the fundamental representation of SU(3). Since 
we do not aim at measuring non-zero momentum glueballs we consider only 
translationally invariant operators, i.e. operators averaged in space. 

It is computationally feasible to measure Wilson loop operators up to length 
8. The composition of the irreducible representations T^^ of the cubic group 
in terms of these 22 loop shapes has been done already in Ref . |Q . Operators 
contribute to the two- and three-dimensional representations with two and three 
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different polarisations, respectively, in analogy to different magnetic quantum 
numbers m for a given angular momentum I in the 0(3) group. Measuring 
all these polarisations may suppress statistical noise more than just increasing 
statistics since the different polarisations of a loop shape are expected to be 
anti-correlated. 

Due to mixing with other states present in pure gauge theory on a periodic 
lattice, such as torelons or glueball pairs, and due to the breaking of the contin- 
uum rotational symmetry by the lattice, the identification and the continuum 
spin assignment of single glueball states are additional vital questions which will 
be addressed in Sections |3.5.5 , 3.5.7 



3.5.3 Simulation Parameters 

To allow for a continuum extrapolation, at least for the lighter states, we decide 
to perform simulations at three different lattice spacings in the range 0.10 fm < 
Os < 0.25 fm in volumes between (1.2 fm)'^ and (2.0 fm)'^. The simulation 
parameters are given in Table The gauge fields are updated by performing 
compound sweeps consisting of 4 pseudo over-relaxation and 1 Metropolis sweep, 
and after every 2 compound sweeps we measure all 22 loop shapes up to length 
8. The measurement is performed on APE-smeared configurations (see Section 
3.2) in order to spatially enlarge the operators, thus improving the overlap 
with the glueball states and reducing unphysical high-momentum fluctuations. 
On the coarse lattice at /3 = 3.15 we use smearing levels Sn,n = 2,4,6,8, 10, 
on the other lattices we use levels iS„, n = 3, 6, 9, 12, 15, always with smearing 
parameter = 0.1. 



f3 




S^xT 


V 


# sweeps / measurements 


3.15 


0.2408(5) 


8^ X 16 


(1.96 fm)3 


241000 / 24100 


3.30 


0.1005(6) 


10^ X 20 


(1.59 fm)3 


99000 / 9900 


3.50 


0.0415(4) 


12^ X 24 


(1.22 fm)3 


115000 / 11500 



Table 12: Run parameters for the glueball simulations using the ^ = 2 perfect 
action. 



3.5.4 Determination of the Energies 

The masses of the lowest states and the first excited states (if possible) in all 
the representations are determined using the variational techniques described 
in Appendix B of p8|. The results in units of the temporal lattice spacing at 
are listed in Tables~p9| - ^ in Appendix ^ together with the number N of 
operators entering into the fitting procedure, the time slices on which the initial 
generalised eigenvalue problem is solved (usually 1/2), the number of operators 
AI in the final fit and the corresponding value of x^/iVoF- These results are 



24 



then multiplied by ■ r^/as from Tables || and @to obtain the glueball masses 
in units of the hadronic scale ro, listed in Table p^. 

In order to obtain reliable estimates of the glueball masses from the varia- 
tional method, one has to pay attention that among the large number of op- 
erators (up to 145 from five smearing levels) there are no correlators entering 
the process that are measured exceptionally bad, i.e. with large errors even at 
small time separation, as this may destabilise the determination. To filter out 
such data, we look at the relative errors of the single operators depending on 
the temporal separation of the creation and annihilation operators and drop all 
correlators whose signal falls below a certain treshold value. Usually, this is 
repeated with different treshold values, yielding different sized sets of operators 
(2 ~ 4 sets per representation). We perform O(IOO) mass determinations on 
each operator set, varying the analysis parameters such as M (see Appendix 
B of 1^) or the fitting range t ~ tmin . . . tmax- It turns out that generally the 
variational method used is stable, provided that the badly measured operators 
are absent from the beginning and provided that no important operators are 
missing. Patently, to find such a window of the number of initial operators, iV, 
may be difficult if there is a rather small total number of operators. 

For some of the heavy states, such as the PC = glueballs, whose masses 

can only be determined from the /3 = 3.50 measurements, we have to resort to 
using the correlators at to = and ti = 1 for the solution of the first generalised 
eigenvalue problem. This has the advantage that the correlators from separa- 
tion to = are positive definite by definition (which is not true for to > 1), 
thus rendering the second generalised eigenvalue problem well defined even if all 
the initial operators are kept, M = N. However, the correlators C{t = 0) are 
under suspicion of containing rather little physical information about the corre- 
lation lengths (which correspond to the decay of the correlators) , actually they 
almost only provide information about the relative normalisation of operators. 
Obviously, operators with a large signal at t = might be preferable because 
their signal/noise ratio is probably better on average, however in general they 
do not correspond to the operators having the largest overlap with the lowest- 
lying states. It turns out that one has to pay much more attention in choosing 
the fit range t — imin ■ • ■ imax when using to = 0, ti = 1 because in this case 
at small t the contamination of the lowest masses due to higher states is much 
more significant. 

During the analysis, it turns out that the mass determinations for the scalar 
representation Af~^ are rather challenging. The ground state receives a larger 
relative error than comparable states where the operators are measured equally 
well; for the first excited state it is very hard to obtain a stable determination, 
and the error of its mass may turn out to be huge (e.g. at /9 = 3.50). Probably, 
it is the underlying vacuum, having the same quantum numbers as the glueball 
state, that is responsible for these troubles. We treat the vacuum just like an- 
other state in this representation, so that the glueball ground state is effectively 
a first excited state and the glueball first excited state is effectively a second 
excited state of the representation. Attempts of using other ways of getting rid 
of the vacuum such as the usual v.e.v. subtraction or the subtraction of large t 
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correlators (that are assumed to contain solely noise, however correlated to the 
noise at lower t) or even more sophisticated methods (like solving a generalised 
eigenvalue problem using large t correlators to dig out the vacuum state) do not 
at all succeed in improving the situation. 



3.5.5 Torelons and Multi-Glueball States 

Besides single glueballs, the spectrum of pure gauge theory on a lattice with 
periodic boundary conditions also contains states consisting of several glueballs, 
torelons or mixed states of glueballs and torelons. Although we expect that 
the operators used to measure glueball energies couple most strongly to the 
single glueball states, other states with similar masses and compatible quantum 
numbers might mix with them and distort the result, i.e. the energy determined 
by the analysis of the correlation matrix may be dominated by a multi-glueball or 
torelon state with smaller energy than the single-glueball state to be measured. 

In principle, there are several means of determining the nature of a state 
that has been measured. Firstly, the simulation may be repeated on lattices of 
different physical size, keeping the lattice spacing fixed. As multi-glueball and 
torelon states show a finite- volume scaling behaviour very different from single 
glueballs, such states stand out and may be dropped from further analysis. 
Secondly, one may measure additional operators that couple strongly to torelon 
or multi-glueball states. Including or excluding these operators in the variational 
method and studying the coefficients obtained from the variational method, the 
mixing strength is determined and states which do not mix considerably with 
any of the additional operators may be safely considered to be single-glueball 
states. 

Finally, using the mass estimates for the low-lying glueballs, one may deter- 
mine the approximate locations of the lowest-lying multi-glueball states. Also, 
the minimal energy of mixing torelons may be estimated using the string for- 
mula, eq. (p3|). 

The first two of the methods mentioned above require additional work and 
computer time and go beyond the scope of this work. The last method, however, 
can be done rather easily. Let us first calculate the energy of the lowest-lying 
torelon states that may interfere with our measurements. Single torelons (see 



Section 3.2) transform non-trivially under Z3 symmetry operations; our opera- 
tors, closed Wilson loops, however, are invariant under these transformations. 
This means that they cannot create single torelon states, however the creation 
of two torelons of opposite center charge is possible. If we assume that they 
do not interact considerably and if our lattice extension is rather large then we 
can use the simple formula E2T ~ 2aL to estimate the energy of a torelon pair 
with momentum zero, wher e a is the string tension (measured by the static qq 



potential, see Section |3.3.l| ) and L is the spatial extent of the lattice. Table 
lists the minimum energies of torelon pairs to be expected on our lattices. 

Note, that a state composed of two opposite center charge torelons and with 
total zero momentum is symmetric under charge conjugation. The operators 
for C = — states therefore do not create such torelon pairs. 
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Channel 


j 


/3 = 3.15 


/? = 3.30 


/3 = 3.50 







2.58(9) 


3.55(13) 


3.65(15) 


A++* 





5.47(53) 


6.83(37) 


6.49(148) 


E++ 


2 


5.63(23) 


5.93(15) 


6.08(28) 


E++* 


2 


7.95(73) 


8.76(48) 


10.66(59) 


T++ 


2 


5.67(25) 


5.82(13) 


6.13(16) 


rr++* 


2 




7.97(69) 


10.16(30) 




3 


6.71(73) 


9.11(64) 


10.92(48) 










14.12(130) 


T++ 


3 


7.54(68) 


9.00(44) 


11.21(33) 


T++* 




8.69(91) 











6.36(59) 


6.13(28) 


6.79(27) 










10.27(110) 


E-+ 


2 


7.41(48) 


8.23(28) 


8.48(29) 


E-+* 




8.35(125) 




12.66(84) 


nn-+ 


2 


7.41(49) 


8.18(21) 


8.57(42) 










15.08(67) 




1 


7.20(39) 


7.69(26) 


9.12(24) 


At- 


3 






11.02(68) 


At-* 








15.55(147) 


rp+- 

^2 




7.75(87) 




10.29(118) 


-^2 








14.98(139) 


E+- 








12.42(42) 


Tr 








10.81(83) 










10.85(83) 


Ar 








12.18(45) 


Ar 








15.34(103) 



Table 13: Final glueball mass estimates in terms of the hadronic scale ro, mcro 
from the measurements using the perfect ^ = 2 action. The continuum spin 
assignment J is given as well. 
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roE2T 



3.15 



10.60 



3.30 



8.24 



3.50 



6.92 



Table 14: The minimum energies of momentum zero torelon pairs on the lattices 
used in the glueball simulations, given in units of r^^ . 



To estimate the lowest energies of multi-glueball states present on our lat- 
tices, we follow the method used by Morningstar and Peardon jl^, described in 
great detail in Appendix F of ||3^ . We assume that the energy of multi-glueball 
states is approximately given by the sum of the energies of the individual glue- 
balls, i.e. that there is no substantial interaction. It is therefore clear that the 
lowest-lying multi-glueball states with zero total momentum are the two-glueball 
states, with energy 



where mi denotes the rest mass of the first glueball with momentum p and 
m2 denotes the rest mass of the second glueball with momentum ~p. The 
masses are taken from our determinations at the respective values of (3 and 
the (lattice) momenta are chosen such that the energy E2G of the glueball pair 
contributing to some representation T^^ is minimised. Note, that glueball pairs 
contribute to different representations of the cubic group, not only depending 
on the representations according to which the single glueballs transform, but 
also on the single glueball momentum p. 

The lower bounds of the multi-glueball energy region as well as the lower 
bound for torelon pairs are indicated in Figures |^-|^, together with the deter- 
mined energies of all the states in the respective representations. It turns out 
that all of the excited states measured and even some of the lowest states in a 
given representation could be affected by torelon pairs or multi-glueball states. 
Of course, it is not at all ruled out that these states are indeed single glueballs, 
however one has to be very careful with the interpretation and has to keep in 
mind that these issues require further study. 

3.5.6 The Continuum Limit 

In order to remove discretisation errors from glueball masses tqhuq obtained at 
finite lattice spacing, we have to perform the continuum limit, i.e. to extrapolate 
the results to = 0. Having obtained results at three values of the coupling 
/3, one of them corresponding to a rather coarse lattice spacing, and having no 
accurate information (e.g. from perturbation theory) about the behaviour of 
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A| A, E T| T, A| A, E T, A, A, E T, A^ A, E T, 
+ + -+ +- 



PC 



Figure 7: The results of the glueball measurements at /9 = 3.15 converted into 
units of ro . The lower bounds of the energies of multi-glueball states (solid lines) 
or torelon pairs (dotted line) are given as well. 



A| A, E T| T, A| A, E T, A, A, E T, A^ 

+ + -+ +- 

PC 



A, E T, T, 



Figure 8: The results of the glueball measurements at P — 3.30 converted into 
units of To . The lower bounds of the energies of multi-glueball states (solid lines) 
or torelon pairs (dotted line) are given as well. 
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A| A, E T| T, A| A, E A, A, E T, A, A, E T, 
+ + -+ +- 



Figure 9: The results of the glueball measurements at /? = 3.50 converted into 
units of To . The lower bounds of the energies of multi-glueball states (solid lines) 
or torelon pairs (dotted line) are given as well. 



the energies depending on the cut-off (about the form of the curve used in the 
extrapolation) this is not an easy task and rather ambiguous. 

We perform two different kinds of continuum limits. Firstly, we extrapolate 
the masses in terms of the hadronic scale Tq, rQiTiQ = meat ■ £, ■ ro/ag, secondly, 
we extrapolate mass ratios racilmc2 of different glueball representations. Gen- 
erally, the procedure we resort to is the following: We include all the three 
energies or mass ratios, obtained aX (i — 3.15, 3.30, 3.50, into the fit and use the 
form 

romc \a, = rorriG |q,=o +C2 f — ) , (38) 

thus including the correction of the smallest power of the (spatial) lattice spac- 
ing Os to be expected to occur. Additionally, we perform fits including only 
the results from the two finer lattices at /3 = 3.30, 3.50 using a constant (if this 
fits the results reasonably well) or the form stated above. This value of the 
continuum mass (or mass ratio) is then compared to the result of the fit of all 
the three data points, which is only accepted if the two results coincide within 
their errors. 
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Masses in Terms of Tq 



For the extrapolation of masses in terms of the hadronic scale ro , this procedure 
works very well for the representations A^"'", E ^, E^~^, Tj""*" and Tj^"*" 

where we may always fit the results from the two finer lattices to a constant 
which agrees with the result of the extrapolation using eq. ( |38| ) on all three data 
points. For the A^^ representation, we obtain the best results from fits with a 
constant to two as well as to three data points. 

Much more difficult is the extrapolation for the scalar glueball and its excited 
state, Af^ and A'^'^ . It turns out that the fits have bad values of and for the 
ground state determinations the result obtained performing a fit to a constant 
mass using the results from the two finer lattices does not correspond to the 
fit including the third, coarser lattice and a quadratic term in a^. We decide 
to keep the (usual) result, including the quadratic term and keep in mind that 
there may be present some ambiguities which will be discussed later. 

The results of the continuum extrapolations for various glueball states are 
collected in Table In Figures ^ and |Tl| the fitting curves are shown together 
with the measured values for the PC = ++ and the PC = — h sectors where 
extrapolations are possible with our data. 

Other representations where energies can be determined at several values 
of (3, such as , T^'^ or T/"'" , do not yield consistent results using the 
procedure presented above and we decide to abstain from performing continuum 
extrapolations in terms of romc for these states. 



Mass Ratios 

Extrapolating glueball mass ratios mci/wca to the continuum has the advan- 
tage that uncertainties present in the determination of ro/ag and the renor- 
malised anisotropy S^r cancel out. Furthermore, we may assume that finite size 
effects are similar for glueball states measured on the same lattice as opposed 
to completely different quantities used e.g. to obtain ro/a^ and ^r. 

The common way of extrapolating glueball mass ratios to the continuum 
is using ratios mQ/mQ++ , i.e. ratios to the scalar glueball mass, which is the 
lightest mass present. Doing this, we observe that the extrapolation is difficult, 
the value of of the fit is rather bad and the errors are large. Furthermore, 
the resulting continuum mass ratios are not in good agreement to ratios of the 
continuum masses extrapolated in terms of tq. This is not a surprise, because 
the behaviour of the scalar glueball as a function of the lattice spacing is rather 
complicated, i.e. there are large cut-off effects that cannot be described too 
easily. Additionally, the errors of the mass estimates of the scalar glueball are 
quite large. Due to these reasons, we decide to use the well measured mass of 
the tensor representation, which show small errors and seem to scale rather 
well. 

The procedure of the extrapolation is the same as in the case of romc de- 
scribed in the previous section; the unambiguous results are listed in Table p3| . 
Again, for the A^'^ states, as well as for the E'^^ , the fit to a constant works 
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out best. Note, that in this way we can obtain rehable continuum results for 
the representations , T^^ and Tj*"^ which is not possible performing the 
extrapolation using masses in terms of tq. Furthermore, all the other results 
listed are in agreement to the ratios of the continuum masses toTOg within one 
standard deviation. 

In the following, we will use the finite lattice spacing results from the simu- 
lation on the finest lattice at /3 = 3.50 for the energies of states that cannot be 
extrapolated to the continuum in one of the two ways, e.g. to study degeneracies 
or to draw a (rather qualitative) picture of the low-lying glueball spectrum. In 
tables listing continuum results, the values not extrapolated to Os = will be 
stated in square brackets. 




Figure 10: Mass estimates of the PC — glueballs in terms of the hadronic 
scale tq against (as/rg)^. The curves are the continuum limit extrapolations 
of the forms as indicated in the text. Circles: ^ boxes: A^^, diamonds: 
E~^^ , upward triangles: T^~^ , downward triangles: T^^; solid symbols: ground 
states, open symbols: first excited states. Note that the continuum results of 
the representations T^^ and T^^ have been obtained using a fit of glueball 
mass ratios. 
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Figure 11: Mass estimates of the PC — — h glueballs in terms of the hadronic 
scale ro against (as/ro)^. The curves are the continuum Umit extrapolations 
of the form romc + C2{cis/i'of' . Circles: A^^ ^ diamonds: E ^, downward 
triangles: T2^] solid symbols: ground states, open symbols: first excited states. 
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3.5.7 Continuum Spin Identification 

Once the extrapolation of the measured glueball energies to the continuum, 
tts = 0, has been performed, what remains is the assignment of continuum spin 



J to the different representations of the cubic group T. From Table 11 we know 
from which lattice representations r^*-^ a continuum glueball with quantum 
numbers J'^'^ may obtain contributions, where the assignment of parity P and 
charge conjugation C is simply one to one. Degeneracies (in the continuum 
limit) between several representations T^*^ contributing to the same continuum 
state J^'~' are a strong indication for the correctness of the assignment of all 
the lattice states involved to the same continuum J^'~^ glueball. Furthermore, 
we make the assumption that the mass of the glueballs increases with the spin 
getting larger. 

In the PC = ++ sector, we observe a single low- lying state, A^'^ , which 
is thus assigned J = 0. The and T^'^ states are degenerate to a very 

high precision as it should happen if these representations correspond to the 
five polarisations of a J = 2 glueball. The excited state A^'^ again has no 
degenerate partner and is assigned to an excited state of the continuum J = 
glueball. The excited states i?"*""*"* and T^~^ again turn out to be degenerate 
and so they are assigned to an excitation of the J — 2 glueball. The same 
situation is met with the representations and T^^ which are thus assigned 
J = 3. Since the first excitation T^'^* is part of the J = 2 state, the missing 
3 polarisations should come from the second excitation of Tj*""*" which, however, 
cannot be measured by our simulations. 

In the PC = — h sector, the situation of the A^'^ and the E ^ and Tj"^ is 
very similar to the one of their partners in the PC — ++ sector; these states are 
thus assigned J = and J — 2 respectively. The remaining three, excited states 
have no apparent degeneracies, there is thus no safe assignment of continuum 
spin to these states. The facts that the excited state A'^'^ is rather light and 
comes with no degenerate partner, raises the presumption that it is the excited 
J = 0"+ glueball. 

In the other two sectors, continuum extrapolations of our measured energies 
are not possible, except for the representation using mass ratios. However, 
looking at the degenerate representations in the PC = ++, — h sectors, it is 
noticeable that the degeneracies are apparent even for finite lattice spacing at 

P = 3.50,3.30. Assuming that this behaviour persists for PC = H — , , at 

least for energies not too high, we may assign continuum spin even to some of 
the remaining states. 

In the PC = H — sector, we notice the lowest-lying state having no 

degenerate partner thence suggesting a J = 1 interpretation. Next, there is the 
A2 which is likely to correspond to J = 3, which is the smallest possible spin 
corresponding to this representation. 

Finally, in the PC = sector, there are different possible scenarios. The 

almost exact degeneracy of T-^ and T2 suggests them contributing to the 
J = 3 continuum state. The rather heavier state could still be degenerate 
with the two latter states (note the high mass and the finite lattice spacing!) 
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and carry the remaining polarisation. Admittedly, it could as well correspond 
to an excited J = 3 state or even to J — 6. The final Al state is very heavy 

and indicates that glueballs with even spin and PC = have large energies 

as the representation Ai contributes solely to even numbered spin states up to 
J = 8. 

Let us remark that a new method of improving the overlap of operators 
to glueballs of a given continuum spin J thus rectifying the continuum spin 
identification is presented in ^ . 

The well supported continuum spin assignments presented above are given 
together with the glueball masses in Tables |l3, 32, 3^, while the final results on 
the masses of these states are collected, together with the masses in MeV, in 
Table ITsl. 







■mcro 


TOG [MeV] 


0++ 


At+ 


4.01(15) 


1580(60) 


0++* 


A++* 


7.66(71) 


3030(280) 


2++ 


E++, T++ 


6.15(15) 


2430(60) 


2++* 


E++*, T++* 


10.52(51) 


4160(200) 


3++ 


At\ T++ 


11.65(49) 


4600(190) 


0-+ 




6.46(18) 


2550(70) 


2-+ 


E-+, Ti+ 


8.75(26) 


3460(110) 


0-+* 




[10.27(110)] 


[4060(430)] 


1+- 


Tr 


9.45(71) 


3730(280) 


3+- 


At- 


[11.02(68)] 


[4350(270)] 



Table 15: Final results for the masses of continuum glueballs with quantum 
numbers J^^ , obtained from the lattice representations F^'-^. For the conversion 
to MeV, ^ sa 0.5 fm w 395 MeV has been used. Quantities in brackets have 
not been extrapolated to the continuum but denote results from the simulation 
on the finest lattice at /3 = 3.50 (corresponding to Us = 0.102 fm). 



3.5.8 Discussion 

In order to discuss the results obtained from the glueball measurements let us 
first list the sources of problems that may affect our results. Firstly, there 
are finite size effects; our finest lattice (/? = 3.50) has got a spatial volume 
of (1.22 fm)"^ which is rather small compared to other volumes used in glueball 
measurements. The magnitude of finite-size effects depends largely on the quan- 
tity studied, more precisely on the form of the glueball wave functions (above 
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all their extension). Due to the three volumes employed having such different 
size, finite-size effects may not only systematically shift the lattice results but 
they can also make reliable continuum extrapolations much more difficult or 
even impossible. 

Secondly, the mass of the scalar glueball A^^ (above all the ground state) 
shows very large cut-off effects. These could be due to the presence of a criti- 
cal end point of a line of phase transitions in the fundamental-adjoint coupling 
plane. Our parametrisation includes in its rich structure operators transform- 
ing in the adjoint representation. If the net coupling of the parametrised action 
(which we do not control during the construction and the parametrisation) lies 
in a certain region, the effect of the critical end-point on scalar quantities at cer- 
tain lattice spacings (sometimes called the "scalar dip" ) may even be enhanced 
compared to other (more standard) discretisations with purely fundamental op- 
erators. It is important to note that the true classically perfect action is not 
expected to be sensitive to this critical point, but it is rather the parametrisation 
which may incidentally cause the sensitivity close to it. 

Furthermore, our analysis shows that states consisting of several glueballs 
or torelon pairs could mix with most of the higher-lying states occurring in 
our measurements. This makes the interpretation of the measured states more 
difficult, additionally, systematic effects (e.g. different mass shifts at different 
values of /? because of the different lattice sizes or due to the presence of op- 
erators having larger overlap with unwanted states) may again complicate the 
continuum extrapolation. 

There are several possibilities of improving the situation. Firstly, and most 
simple, improving the statistics may help. This can be seen e.g. by looking at 
the results of the simulation on the coarsest lattice at /9 = 3.15. Although the 
lattice is coarser, more energies can be determined than on the finer f3 — 3.3- 
lattice, mainly due to the larger statistics. Additionally, a lot of problems stated 
previously could be kept under control if for each value of (3, we determine the 
glueball energies on several lattices of different physical size. A way of improving 
the measurements is systematically studying the operators employed to create 
and annihilate single glueball states in order to increase their overlap with the 
single glueballs to be measured and to decrease their overlap with all the other 
states. This includes the study of other smearing techniques, such as e.g. Teper 
fuzzing . Conversely, the introduction of operators coupling most strongly 
to the unwanted multi-glueball or torelon-pair states allows for the calculation 
of mixing strengths in order to exclude all the unwanted states from further 
analysis. Another technically rather easy but computationally expensive way of 
improving the results is performing simulations on additional values of /3. This 
allows for a more detailed study of the continuum limit, especially in the case 
of non-perturbatively improved actions, where there is no clearcut information 
about the way how the continuum is approached. Naturally, one has to take 
care not to go too deep into the strong coupling region as there is not much 
information about the continuum. 

There is one error which is not cured by the measures recommended in the 
previous paragraph, namely the effects on the scalar states, coming from the 
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"scalar dip" . Morningstar and Peardon improve the situation by adding a two- 
plaquette adjoint term with a negative coefficient which results in an approach 
to the continuum on a trajectory always far away from the dangerous "dip" 
region Q. In principle, the classically perfect action could be treated the 
same way: extract all the operators with adjoint contributions present in the 
parametrisation, determine their sign and add another (non-linear) constraint 
to the fit, namely that the action of all these operators together corresponds 
to an adjoint operator with a certain (negative) coupling. Because of the large 
freedom in the fit, the inclusion of this single criterion should not impair the 
quality of the parametrisation considerably. 



Collab. 




r 01712+ + 


UKQCD I?) 


4.05(16) 


5.84(18) 


Teper @ 


4.35(11) 


6.18(21) 


GFll |6| 


4.33(10) 


6.04(18) 


M&P 


4.21(15) 


5.85(8) 


Liu [H 


4.23(22) 


5.85(23) 


FP action 


4.12(21) 


[5.96(24)] 


AICP 


4.01(15) 


6.15(15) 



Table 16: Comparison of the lowest-lying glueball masses in units of tq. Values 
in brackets denote masses obtained at a lattice spacing a = 0.10 fm and are not 
extrapolated to the continuum. 



Collab. 



romQ-+ 



romo++* 



rom2-+ 



romi+- 



Teper | 
M&P Il3 



5.94(68) 
6.33(13) 



7.86(50) 
6.50(51) 



8.42(78) 
7.55(11) 



7.84(62) 
7.18(11) 



FP action 
AICP 



[6.74(42)] 
6.46(18) 



7.66(71) 



[8.00(35)] 
8.75(26) 



[7.93(78)] 
9.45(71) 



Table 17: Comparison of glueball masses in units of vq. Values in brackets denote 
masses obtained at a lattice spacing a = 0.10 fm and are not extrapolated to 
the continuum. 



Table |l6| and |l7| compares our continuum glueball masses obtained with the 
anisotropic classically perfect action (AICP) to the results of the isotropic FP 
action as well as to results obtained by other collaborations. There is reasonable 
agreement between the different determinations of the 0++, 2+"'", 0+"'"*, '", 
2 glueball masses. The situation for the 1^' state is different, however. This 
is the heaviest state that occurs in our analysis and correspondingly difficult to 
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determine on the lattices used, having a rather small anisotropy ^ = 2. Thinking 
of all the possible sources of errors stated above, which are (partly) also present 
in the analyses of the other groups, we tend to explain this discrepancy with 
underestimated (or disregarded) systematic errors as discussed previously. In 
particular, by looking at the finite lattice spacing results for in Table (p^, 
we observe that the masses determined on the coarser (and larger) lattices at /? = 
3.15, 3.30 coincide with the (continuum) results obtained by other groups and 
do not even show a significant discrepancy, despite the considerable difference 
between the lattice spacings. The mass determined on the fine (and small) 
lattice at /3 = 3.50 however, is much higher. This tendency is even amplified 
by the continuum extrapolation. We suspect that this state exhibits strong 
finite-size effects pushing up its mass in small volumes. 




Figure 12: Lattice results and continuum extrapolation of the scalar (O'^''^) 
glueball mass for the anisotropic classically perfect action, together with results 
obtained with different other actions. The solid and dotted lines are the con- 
tinuum extrapolations of the AICP and M&P data, respectively. The different 
continuum values for the Wilson action stem from different groups. 

Figures ^ and ^ compare our measurements of the , and Tj^^ 

states to measurements obtained by other groups, using different actions, as 
well as to the isotropic FP action results. Concerning the scalar glueball, the 
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Figure 13: Lattice results and continuum extrapolation of the tensor (2++) 
glueball mass for the anisotropic classically perfect action, together with results 
obtained from simulations employing different other actions. The continuum 
extrapolations of the E^'^ and T^~^ representations are averaged in order to 
get the mass of the single continuum glueball. The solid and dotted lines are 
the continuum extrapolations of the AICP and M&P data, respectively. The 
different continuum values for the Wilson action stem from different groups. 
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artifacts of the isotropic FP action at moderate lattice spacings corresponding 
to (os/ro)^ < 0.15 {us < 0.19 fm) seem to be rather smaher than the ones of the 
tadpole and tree-level improved action (M&P) and certainly much smaller than 
the artifacts of the Wilson action. For the AICP action, however, the situation is 
less clear. The mass obtained on the coarsest lattice at /3 = 3.15 (corresponding 
to Us « 0.24 fm) exhibits large cut-off effects of about 35% (compared to about 
20% for the M&P action), probably due to the scalar dip. Concerning the tensor 
glueball, the picture is not clear, mainly due to the considerable statistical errors 
of the measurements, except for the ones obtained using the tadpole and tree- 
level improved anisotropic action. 



4 Repeating the Spatial Blocking Step 

As noted in section |[ the spatial blocking step used to obtain a classically 
perfect ^ = 2 action may be repeated straightforwardly in order to generate 
actions with higher anisotropics. Using another spatial scale 2 blocking step we 
create a classically perfect C = 4 gauge action and parametrise it in the same way 
as the ^ = 2 action. Furthermore, we measure the renormalised anisotropy and 
show that it is indeed feasible to construct actions suitable for MC simulations 
of heavy states such as glueballs. The classically perfect ^ = 4 action is not yet 



examined thoroughly, however the parametrisation is ready (see Appendix A. 5) 
and the analyses of the ^ = 2 action (see section ||) may be repeated for ^ = 4. 

4.1 Construction 

Each application of the spatial blocking step (slightly) renormalises the anisotro- 
py of the action; the renormalised anisotropics of the final actions thus have to 
be measured in the end. However, this is necessary anyhow if one is interested in 
comparing most results (other than e.g. mass ratios, where the renormalisation 
of the anisotropy cancels out) to other collaborations or to the experiment. 

In order to be able to repeat the spatial blocking step getting from ^ = 2 
to ^ = 4 we construct a ^ = 2 action which is valid on minimised ^ = 4 con- 
figurations including into the fit configurations at different values of /3 over a 
large range. The ^ = 2 action presented in section |^ is not suitable for this task 
as it is dedicated to be used solely on largely fluctuating configurations around 
P = 3.0. The construction and the parameters of a suitable intermediate ^ = 2 



action are described in Appendix A. 4 



Once the fine action to be used on the r.h.s. of the renormalisation group 
equation, eq. (0), is ready, we may proceed along the same lines as for the 
^ = 2 action. We perform the non-linear fit on 2 ^ad-hoc ~ 6 configurations 
each at /3ad-hoc = 4.0, 3.5, 3.0, 2.5, 2.0. It turns out that in order to obtain 
parametrisations that are free of dangerous "traps" in the u-w plane (see Section 



3.1) we have to include the condition A{u,'w) > at ten points {u,w) into the 



non- linear fit. After studying the values of a-s well as the linear behaviour 
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of the parametrisations we decide to use a set with inax(fc + /)sp = 3 and 
inax(fc + Z)tm = 2. The parameters are given in Appendix A. 5. 



4.2 The Renormalised Anisotropy 

To check whether the construction of the classicahy perfect ^ = 4 action works 
well and whether the paramctrisation reproduces the input anisotropy, we mea- 
sure the renormalised anisotropy at one value of /3 = 3.0, using the torelon 



dispersion relation, following the method described in Section 3.2. The simula- 
tion parameters are collected in Table 





S"^ X LxT 


# sweeps / measurements 


3.0 


8^ X 4 X 32 


54000 / 10800 



Table 18: Run parameters for the torelon measurements using the ^ = 4 perfect 
action. The lattice extension in torelon direction L, the extension in the two 
transversal spatial directions S as well as the temporal extension L are given. 



Due to the (expected) coarser spatial lattice spacing we do not need as many 
smearing steps as for the = 2 action and thus perform our measurements on 
smearing levels 5„,n = 2, 4, 6, 8, 10 keeping =0.1 fixed. 

Figure |l^ displays the dispersion relation so obtained, while Table ^ col- 
lects the measured energies E{p'^) determined using variational methods (see 
Appendix B of |p8|), together with the number of operators used in the vari- 
ational method, the fit ranges and the values of P^r degree of freedom, 
X^/A'df- The renormalised anisotropy as well as the torelon mass are deter- 
mined in exactly the same way as for the ^ = 2 action. The results are listed in 
Table |l^. Furthermore, we may again evaluate an estimate of the lattice scale 
using the string picture relation, eq. (|30|), yielding — 0.532(16), and thus 

fo/cis = 2.241(86) which corresponds to = 0.22 fm, at — 0.060 fm. Having 
a torelon of length Lag « 0.89 fm the error of the scale estimate is expected to 
be about 5-10% (see Table |). 



/3 


fit range 




TUT at 


xVNuf 


3.0 


1..5 


3.71(8) 


0.235(13) 


0.82 



Table 19: Results of the torelon simulations using the ^ = 4 perfect action. The 
fit range in p"^ is given in units of {2tt/S)^. 



Concerning the repeated application of the spatial blocking method pre- 
sented in this work to yield classically perfect gauge actions for higher anisotropics, 
we may conclude that the construction including the blocking as well as the 
paramctrisation works in exactly the same way as the construction of the £, — 2 
action. The resulting action shows a renormalisation of the anisotropy of about 
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1 q 




1 1 < ' < ' < ' < L 

2 4 6 8 

2 2 

(a p) in units of (2ii/8) 



Figure 14: Torelon dispersion relation for (3 = 3.0. The straight Hne is the 
correlated fit to E'^{p) = m% in the range = 1..5. 
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7% at the value of /? = 3.00 (corresponding roughly to ag ~ 0.2 fm). Calcu- 
lations of physical observables such as the static quark-antiquark potential or 
glueball masses using the ^ = 4 classically perfect action are presently absent. 

5 Conclusions 

In this work, we have presented the construction and parametrisation of a clas- 
sically perfect anisotropic SU(3) gauge action based on the Fixed Point Action 
technique. The recently parametrised isotropic FP action and its parametri- 
sation ansatz using mixed polynomials of plaquettes built from single gauge 
links as well as from smeared links have been the starting point. Performing 
one and two purely spatial blocking (RG) steps, respectively, starting from the 
isotropic action, and adapting the parametrisation appropriately, we have ob- 
tained parametrised classically perfect actions with anisotropics ^ = 2 and 4. 
The ^ = 2 action has been tested extensively in measurements of the torelon 
dispersion relation, of the static quark-antiquark potential, of the deconfining 
phase transition as well as of the low-lying glueball spectrum of pure gauge 
theory. 

The results of the torelon measurements show that the renormalisation of 

the anisotropy (due to quantum corrections and parametrisation artifacts) is 
small (below 10% for ^ = 2 and 4). The rotational invariance of the action has 
been examined by measuring the static quark-antiquark potential, including 
separations corresponding to a large set of lattice vectors, again with reassuring 
results. The glueball measurements, including some rather heavy states, confirm 
that the use of anisotropic lattices facilitates spectroscopy when heavy states are 
present. Compared to the isotropic simulations, with the same amount of com- 
putational work it is possible to resolve states with considerably larger energy, 
allowing continuum extrapolations from larger ranges of the lattice spacing. 
The mass of the scalar ghieball with quantum numbers J^'^ = is mea- 
sured to be comparable (1580(60) MeV) to masses obtained by other groups 
(around 1670 MeV); however it shows large cut-off effects when measured on 
a lattice with spatial lattice spacing ~ 0.25 fm. This could be caused by 
a large sensitivity of the action to a critical point in the fundamental-adjoint 
coupling plane, due to adjoint terms in the action. Additionally, one has to be 
aware of the manifold sources of possible errors in measurements of the glueball 
spectrum. Having at hand the statistics reached so far, no ultimate statement 
about the goodness of the parametrised classically perfect actions can be made 
yet. 

During our studies of both isotropic and anisotropic classically perfect ac- 
tions, we have noticed that the examination of scaling properties of lattice gauge 
actions is a very delicate problem: Quantities which can be reliably measured 
and which are not very sensitive to systematic factors (like the volume, other 
states mixing with the observed ones or the method used for extracting the 
mass etc.), e.g. the critical temperature or the hadronic scale rg (for moderate 
lattice spacings) generally exhibit rather small cut-off effects and thus demand 
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very large statistics, if the differences between the actions should be explicitly 
identified. On the other hand, some glueball states show large cut-off effects 
which makes them interesting objects to study scaling violations. As pointed 
out in this work, however, there are a lot of systematic factors making the ex- 
traction of lattice artifacts difficult. The task of comparing different actions is 
even more difficult if the measurements and analyses are performed indepen- 
dently by different groups as this introduces additional systematic discrepancies 
and ambiguities. 

Let us therefore conclude that despite the rich parametrisation which allows 
for all the beautiful properties of classically perfect actions like scale invariance 
of instanton solutions or rotational invariance, there is no conclusive evidence 
of the parametrised classically perfect actions (isotropic as well as anisotropic) 
behaving significantly better than other improved actions (such as the tadpole 
and tree-level improved anisotropic gauge action) — however, there is also no 
evidence for the converse. Certainly, one has to consider the large overhead of 
the parametrisation as compared to other actions if one is about to plan pure 
gauge simulations. 

For the future, very accurate scaling tests comparing the classically perfect 
actions to other improved gauge actions are desirable. On one hand, these 
can include large statistics measurements of the critical temperature and the 
hadronic scale, on the other hand, one might perform extensive simulations of 
the glueball spectrum, including all the limits to be taken and systematically 
excluding all known sources of errors. 

With respect to the glueball measurements another promising plan would 
be to gain control about the adjoint operators present in the parametrisation 
of the classically perfect actions in order to circumvent the critical point in the 
adjoint coupling plane such that the influence of the scalar dip is minimised, 
similarly to what has been done for the tree-level and tadpole improved actions 
by Morningstar and Peardon |Q. 
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A Action Parameters 

In this Appendix we collect the parameters of the FP and the classically perfect 
actions that have been used or constructed throughout this work. The isotropic 
FP action for coarse configurations has been presented already in pS|. The 



isotropic action in Section A.l is an intermediate parametrisation valid on con- 
figurations typical for MC simulations that are minimised once, parametrised 
during the cascade process for the isotropic action. It is used in the spatial block- 



ing (see Section pAh to minimise the coarse S, — 2 configurations, constructed 



using the ad-hoc anisotropic action presented in Section A. 2. The resulting 



^ = 2 perfect action for coarse configurations is presented in Section A. 3 
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To repeat the blocking step, we need a ^ = 2 action which is valid for ^ = 4 
configurations minimised once in a purely spatial blocking step. This action is 
described in Section A. 4. The resulting ^ = 4 action, finally, is presented in 



Section A. 5 



A.l The Intermediate Isotropic Action 

The intermediate isotropic action has been parametrised during the construc- 
tion of the isotropic FP action. It is supposed to be valid on configurations that 
are obtained by minimising configurations typical for MC simulations. It fulfills 
the O(a^) Symanzik conditions (see |2^) and uses polynomials in the fluctu- 
ation parameter x^in), thence it smoothly approaches the continuum limit and 
is expected to interpolate between the rather coarse configurations mentioned 
above and the smooth limit. It is not intended to be used in MC simulations. 



since its linear behaviour (see Section 3.1) is not checked 



This action has been used for the (spatial) minimisation of coarse ^ = 2 
configurations, i.e. to describe the isotropic configurations on the r.h.s. of eq. (^. 

The non-linear parameters describing polynomials of order 3 are (cf. formula 
(i, § and (0), (H): 



^(0) 


= 0.082, 


^(1) 


= 0.292353, 




= 0.115237, 


^(3) 


= 0.011456, 


cT 


= 0.282, 




= -0.302295, 


^1 


= -0.302079, 


''1 


= -0.052309, 




= 0.054, 


^2 


= 0.298882, 


J2) 
C2 


= -0.081365, 


Si) 
C2 


= -0.023762, 


4°^ 


= -0.201671, 


^3 


= 0.022406, 




= 0.004090, 




= 0.014886, 


4°' 


= -0.008977, 


„(1) 
C4 


= 0.245363, 


C4 


= 0.140016, 


C4 


= 0.028783. 



The linear parameters are collected in Table 



Pkl 


I = 


I = 1 


/ = 2 


I = 3 


I = 4 


k = 







0.629227 


-0.556304 


0.186662 


-0.010110 


k = 


1 


-0.368095 


0.852428 


-0.199034 


0.031614 




k = 


2 


0.389292 


-0.207378 


-0.010898 






k = 


3 


-0.054912 


0.039059 








k = 


4 


-0.000424 











Table 20: The linear parameters pki of the parametrised intermediate isotropic 
FP action. 



A. 2 "Ad-hoc" Anisotropic Actions 

In the spatial blocking procedure described in section]^ one needs coarse anisotropic 
(^ = 2, 4, 6, . . . ) gauge configurations which are spatially minimised leading to 
^' — £,/2 configurations. As the goal of this step is to obtain a perfect anisotropic 
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action with anisotropy ^ these coarse anisotropic configurations have to be pro- 
duced using some other action that is ah'eady present. This requirement might 
seem to endanger the whole ansatz, however it is not crucial how the coarse 
configurations exactly look like as the perfectness of the resulting coarse action 
comes from the perfect action on the fine configuration as well as from the exact- 
ness of the RG transformation. Still, we try to create coarse configurations that 
might look similar to future ensembles produced using the perfect anisotropic 
action on the coarse level and whose minimised configurations appear to have 
an anisotropy approximately ^' — ^/2. 

In order to achieve this, we modify the isotropic FP action by adding a term 

— (where pf q denotes the simple temporal plaquette) . This modification 
turns the isotropic Wilson action into the Wilson action with bare anisotropy 
^ and is expected to work approximately also for our FP action. The main 
argument of using this coarse action and not e.g. the anisotropic Wilson action 
is that due to the spatial lattice spacing as being larger than the temporal 
spacing at, the O(a^) artifacts are also expected to be larger than the 0(0^) 
effects. The modification described above should (approximately) preserve the 
FP properties corresponding to Us and is thus expected to be considerably better 
than the naive anisotropic Wilson action. 

Using this ad-hoc modification of the isotropic FP action ( "ad- hoc" action) 
in MC simulations shows (as expected) that (at least for small anisotropies 
^ < 5) its properties on the glueball spectrum are not comparable with the ones 
of the isotropic FP action, however the generated ensembles resemble to the 
ones generated with true perfect anisotropic actions much more than ensembles 
generated with the Wilson action. 

The anisotropies ^ad-hoc that have to be used to generate coarse anisotropic 
configurations turning into minimised configurations with anisotropy ^' are 
^ad-hoc ~ 3.2 for ^' = 1 and ^ad-hoc ~ 6 for ^' = 2. However, this value of 
^ad-hoc varies considerably with /3 — but as stated at the beginning of this 
section, the exact form of the coarse configurations is not essential. 

A.3 The ^ = 2 Perfect Action 

The ^ = 2 perfect action uses the parametrisation described in Section ^j^. The 
number of non-zero asymmetry values rif'^ is 4, the parameters cf^^ (? = !,... ,3) 
are splitted into 3 parameters depending on the contribution to the smeared 
plaquette. The linear parameters pki are non-zero for < fc -I- Z < 4 for spatial 
plaquettes and < A: -I- ^ < 3 for temporal plaquettes. 

The non-linear parameters (constants in x^) have the values 

r?i = -0.866007, r?2 = -0.884110, 773 = 2.212499, 774 = 1.141177; 

cii = 0.399669, C12 = 0.519037, C13 = -0.071334, 

C21 = -0.076357, C22 = -0.031051, C23 = -0.282800, 

C31 = 0.032396, C32 = -0.015844, C33 = -0.046302. 
The linear parameters are collected in Table ^ 
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pU 


7 r\ 
I = U 


7 -1 
1 = 1 


7 O 
1 = 1 


7 Q 
t = O 


7 A 
/ = 4 


k = 







0.433417 


0.098921 


-0.116251 


0.023295 


k = 


1 


0.217599 


-0.272668 


0.248188 


-0.045278 




k = 


2 


0.316145 


-0.180982 


0.028817 






k = 


3 


-0.039521 


0.003858 








k = 


4 


0.005443 










Pl\ 


; = o 


I = 1 


I = 2 


I = 3 




k = 







-0.190195 


0.554426 


-0.121766 




k = 


1 


1.521212 


-0.328305 


0.086655 






k = 


2 


0.011178 


0.020932 








k = 


3 


0.022856 











Table 21: The linear parameters of the ^ = 2 parametrised classically perfect 
action. 

A. 4 The = 2 Intermediate Action 

In order to be able to repeat the spatial blocking step constructing a ^ « 4 
action based on the ^ = 2 perfect action we need a parametrisation of the ^ = 2 
action which is valid on » 2) configurations that arc obtained by spatially 
minimising coarse ^ = 4 configurations once. To construct such an action, we 
perform a non-linear fit to the derivatives of 5 sets of two configurations each 
at /3ad-hoc = 6, 10, 20, 50, 100. The non-linear parameters arc chosen to be 
linear in the fluctuation parameter Xij,{n). Having four different parameters 
T] and splitting up Ci into three parameters (as it is done for all anisotropic 
paramctrisations), this makes 20 non- linear parameters to be fitted which is 
quite at the edge of what is still possible on our computers and that is why we 
restrict the total number of configurations to 10. Rough checks performed on 
a larger number of configurations, with an even larger number of parameters 
show however, that the resulting non-linear parameters are stable and describe 
the data accurately. 

The derivatives and action values of 5 sets of 10 configurations each are 
included in the linear fit (where the relative weight of the action values is chosen 
to be 1.9 • 10~^ for the configurations at /3ad-hoc = 50, 100 and 7.6 • 10~^ at /? = 
6, 10, 20). A linear set where the parameters pki are non-zero for < A; + Z < 3 
for spatial plaquettes and < A; -|- / < 4 for temporal plaquettes describes the 
full action very well concerning this data. Again, this parametrisation is not 
intended to be used in MC simulations, thus the linear behaviour of the action 
is not checked. 

The non- linear parameters (up to first order in x^) have the values 
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r?f ^ = -1.861267, 
?7^°> = 2.750293, 
4°' = 0.520960, 
c^a' = 0.266240, 

Cl3 



0.159372, 





= -0.327466, 




-1.075610, 




= -0.550398, 




= 0.089874, 




1.107017, 




= 0.265817; 


•^11 


= 0.006339, 


^21 ^ 


-0.075219, 


Jl) 

'-21 


= 0.059506 


'-12 


= 0.121035, 


^22 ^ 


-0.080771, 


^-22 


= -0.021515 


„(1) 
^13 


= 0.039564, 


''23 — 


-0.043901, 


'-23 


= 0.009672 



The linear parameters are collected in Table 





1 = 


/ = 1 


I = 2 


/ = 3 




k = 







0.088016 


0.002225 


-0.000285 




k = 


1 


0.341850 


-0.015888 


-0.004087 






k = 


2 


-0.053007 


0.010121 








k = 


3 


0.010500 












I = 


I ^ 1 


/ = 2 


/ = 3 


I = 4 


k = 







0.280043 


5.077727 


-13.714872 


12.739964 


k = 


1 


1.343946 


-6.934825 


27.673937 


-32.288928 




k = 


2 


2.069084 


-17.392027 


28.248910 






k = 


3 


3.691733 


-9.584760 








k = 


4 


0.712244 











Table 22: The linear parameters of the intermediate ^ = 2 parametrised classi- 
cally perfect action. 



A. 5 The = 4 Perfect Action 

The ^ = 4 perfect action uses the parametrisation described in Section The 
number of non-zero asymmetry values vif^^ is 4, the parameters c^*^^ (i = 1, . . . ,3) 
are splitted into 3 parameters depending on the contribution to the smeared 
plaquette. The linear parameters pki are non-zero for < fc + Z < 3 for spatial 
plaquettes and < A: -I- / < 2 for temporal plaquettes. 

The non-linear parameters (constants in Xp) have the values 

771 = -1.491457, 772 = -1.115141, riz = 1.510985, ?74 = 7.721347; 

cii = 2.014408, C12 = 0.128768, C13 = 0.162296, 

C21 = -0.915620, C22 = 0.134445, C23 = -0.013383, 

C31 = 1.166289, C32 = 0.061278, C33 = 0.000759, 
The linear parameters are collected in Table 
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Pit 


1 = 


/ = 1 


1 = 2 


/ = 3 


k = 







0.027625 


0.000052 


0.000000 


k = 


1 


0.072131 


-0.016852 


-0.000054 




k = 


2 


0.036818 


0.003558 






k = 


3 


-0.007413 








Pkl 


/ = 


1 = 1 


1 = 2 




k = 







0.795779 


0.621286 




k = 


1 


2.130563 


-0.286602 






k = 


2 


0.076086 









Table 23: The linear parameters of the ^ = 4 parametrised classically perfect 
action. 
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B Collection of Results 







M 


fit range 






3.0 





2 


1-5 


1.368(30) 


0.19 




1 


3 


1-6 


1.372(17) 


0.49 




2 


3 


1-6 


1.451(17) 


0.28 




4 


3 


1-8 


1.570(21) 


0.46 




5 


3 


1-6 


1.610(19) 


0.34 




8 


2 


1-3 


1.707(29) 


0.05 


3.15° 





3 


2-5 


0.688(16) 


0.11 




1 


4 


1-5 


0.806(5) 


0.06 




2 


4 


1-5 


0.902(6) 


0.59 




4 


3 


1-5 


1.059(8) 


2.35 




5 


4 


1-5 


1.151(7) 


0.16 




8 


4 


1-5 


1.310(15) 


0.71 




9 


4 


1-5 


1.356(23) 


1.58 




10 


4 


1-3 


1.470(17) 


0.42 




13 


3 


1-3 


1.621(23) 


0.01 




18 


3 


1-3 


1.762(47) 


0.56 


3.15'' 


1 


3 


1-5 


1.418(23) 


0.05 




2 


3 


1-5 


1.541(37) 


0.53 




4 


3 


1-4 


1.751(67) 


1.20 




5 


3 


1-5 


1.922(104) 


1.08 



Table 24: Collection of results of the torclon measurements using the ^ = 2 
perfect action. For each /3- value and momentum ~ +Py we list the number 
of operators M kept after the first truncation in the variational method, the 
plateau region on which the fit of the correlators to the form Z{p^) exp{—tE{p'^)) 
is performed (fit range), as well as the extracted energy together with the 
per degree of freedom (x^/-/Vdf)- {Continued in Table |i|) 
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3.3 
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1-10 


0.318(6) 
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is 


4 


i— b 
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r 




i-D 
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1 /If; 
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4 


3-7 


A oo 1 / r \ 

0.321(5) 


1 AA 

1.00 


2 


5 


1-8 


0.432(3) 


0.81 


4 


5 


1-8 


0.586(4) 


0.29 


5 


5 


1-7 


0.656(3) 


1.24 


8 


5 


2-7 


0.806(14) 


0.36 


9 


5 


1-7 


0.859(7) 


0.74 


10 


5 


1-8 


0.913(4) 


0.68 


13 


5 


2-5 


1.010(16) 


0.32 


18 


3 


1 5 


1.215(12) 


0.17 



Table 25: Collection of results of the torelon measurements using the ^ = 2 
perfeet action. For each /3- value and momentum — we list the number 

of operators M kept after the first truncation in the variational method, the 
plateau region on which the fit of the correlators to the form ZijP') exp{—tE{p^)) 
is performed (fit range), as well as the extracted energy together with the 
per degree of freedom (x^/-/Vdf)- {Continuation from Table 
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r 


\t\ 


M 


fit range 


atV{r) 


XV^DF 


(1,0,0) 


1 


3 


3-10 


0.3006(1) 


0.20 


(1,1,0) 


1.414 


3 


3-10 


0.3814(2) 


0.21 


(1,1,1) 


1.732 


3 


3-10 


0.4250(3) 


0.65 


(2,0,0) 


2 


3 


6-10 


0.4469(6) 


0.22 


(2,1,0) 


2.236 


3 


4-10 


0.4752(4) 


0.41 


(2,1,1) 


2.449 


3 


4-10 


0.4971(5) 


0.83 


(2,2,0) 


2.828 


3 


4-10 


0.5309(7) 


0.59 


(3,0,0) 


3 


3 


4-10 


0.5445(7) 


1.71 


(2,2,1) 


3 


3 


5-10 


0.5456(8) 


0.62 


(2,2,2) 


3.464 


3 


4-10 


0.5852(9) 


1.53 


(4,0,0) 


4 


3 


4-10 


0.6269(1) 


0.51 


(3,3,0) 


4.243 


3 


5-10 


0.6455(16) 


0.94 


(4,2,0) 


4.472 


3 


5-10 


0.6630(14) 


0.86 


(4,2,2) 


4.899 


3 


5-10 


0.6964(18) 


0.34 


(5,0,0) 


5 


2 


5-10 


0.7044(22) 


1.24 


(3,3,3) 


5.196 


3 


4-10 


0.7215(17) 


1.31 


(4,4,0) 


5.657 


2 


6-10 


0.7509(41) 


1.23 


(4,4,2) 


6 


3 


5-10 


0.7783(25) 


0.77 



Table 26: Collection of results of the off-axis qq measurements using the ^ = 2 
perfect action at /? = 3.30. For each separation vector we list the length of the 
vector in spatial lattice units, the number of operators M kept after the first 
truncation in the variational method, the plateau region on which the fit of the 
correlators to the form Z{r) exp{—tE{r)) is performed (fit range), as well as the 
extracted energy together with the per degree of freedom (x^/-^df)- 
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r 


\A 


M 


fit range 


atV(f) 




(1,0,0) 


1 


3 


2-7 


0.5053(3) 


0.70 


(1,1,0) 


1.414 


3 


3-5 


0.7322(7) 


0.01 


(1,1,1) 


1.732 


3 


3-8 


0.8902(14) 


0.60 


(2,0,0) 


2 


3 


2-8 


0.9407(10) 


0.89 


(2,1,0) 


2.236 


3 


2-8 


1.0642(11) 


0.36 


(2,1,1) 


2.449 


2 


2 5 


1.1676(13) 


0.27 


(2,2,0) 


2.828 


3 


2-6 


1.3082(18) 


0.17 


(3,0,0) 


3 


3 


2-6 


1.3392(25) 


0.16 


(2,2,1) 


3 


3 


2-5 


1.3867(20) 


0.26 


(2,2,2) 


3.464 


3 


2-6 


1.5675(55) 


0.84 


(4,0,0) 


4 


2 


1-5 


1.7429(24) 


0.14 


(4,2,0) 


4.472 


3 


2-6 


1.9328(93) 


1.47 


(4,2,2) 


4.899 


3 


1-4 


2.1197(34) 


0.15 



Table 27: Collcetion of results of the off-axis qq measurements using the ^ = 2 
perfect action at /3 = 3.00. For each separation vector we list the length of the 
vector in spatial lattice units, the number of operators M kept after the first 
truncation in the variational method, the plateau region on which the fit of the 
correlators to the form Z{f) exip{—tE{r}) is performed (fit range), as well as the 
extracted energy together with the per degree of freedom (x^/-/Vdf)- 
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r 


M 


fit range 


atV{r) 


XV^DF 


3.15 


1 


5 


3 10 


0.3703(2) 


0.456 




2 


5 


3-7 


0.6110(4) 


0.735 




3 


5 


3-10 


0.8019(8) 


1.209 




4 


5 


3-7 


0.9775(16) 


1.238 




5 


5 


4-9 


1.129(8) 


0.519 


3.50 


1 


5 


3-12 


0.25910(5) 


1.163 




2 


5 


3-10 


0.36490(12) 


0.769 




3 


5 


7 12 


0.42168(38) 


0.664 




4 


5 


5-11 


0.46587(51) 


0.205 




5 


5 


5-12 


0.50388(76) 


1.139 




6 


5 


5-11 


0.54051(96) 


0.853 



Table 28: Collection of results of the on- axis qq rneasnrerncints using the ^ = 2 
perfect action at /3 = 3.15, 3.50. For each separation along the axes we list 
the number of operators M kept after the first truncation in the variational 
method, the plateau region on which the fit of the correlators to the form 
Z{r)ex.p{—tE{r)) is performed (fit range), as well as the extracted energy to- 
gether with the P^r degree of freedom (x^/iVop)- 
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Channel N 



to/ti M tit range 



XV^DF 



energies 



A++ 
^1 


O i 




1 


1 /I 


U.Dy 




A++ 


59 


1/2 


10 


1-3 


1.23 


1.365(104) 


E++ 


47 


1/2 


7 


1 - 3 


0.10 


1.405(29) 


E++* 


47 


1/2 


5 


1 - 3 


0.17 


1.985(140) 


rp + + 


47 


1/2 


10 


1 - 4 


0.10 


1.416(32) 


rr++* 


22 


1/2 


6 


1 - 3 


1.20 


1.990(131) 




21 


1/2 


2 


1 - 3 


0.03 


1.675(147) 


T++ 


36 


1/2 


4 


1 ~ 3 


0.34 


1.881(131) 


T++* 


36 


1/2 


6 


1 3 


0.31 


2.170(182) 


AI^ 


25 


1/2 


3 


1 - 3 


0.001 


1.588(114) 


E-+ 


25 


1/2 


2 


1 3 


0.01 


1.849(81) 


E-+* 


25 


1/2 


2 


1 - 3 


2.19 


2.084(268) 




64 


1/2 


4 


1 - 3 


0.02 


1.850(83) 




61 


1/2 


4 


1 - 3 


0.24 


1.798(59) 




18 


1/2 


4 


1 - 3 


0.41 


1.935(177) 



Table 29: Results from fits to the ^ = 2, /3 = 3.15 glueball correlators on the 
8^ X 16 lattice in units of the temporal lattice spacing at- to/ti are used in the 
generalised eigenvalue problem, N is the number of initial operators measured 
and M denotes the number of operators kept after the truncation in C{to)- 
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Channel 


to/ti 


M 


fit range 




energies 


Ar 


1/2 


9 


2-7 


0.90 


0.590(17) 


A++* 


1/2 


11 


1 - 3 


1.57 


1.133(53) 


E++ 


1/2 


11 


1 - 5 


0.97 


0.983(17) 


E++* 


1/2 


8 


1 - 3 


0.90 


1.453(68) 


T++ 


1/2 


12 


1 -4 


1.09 


0.965(15) 


r++* 


1/2 


10 


1 -4 


1.68 


1.386(51) 




1/2 


4 


1 - 3 


0.01 


1.511(94) 


T++ 


1/2 


8 


1 - 5 


0.84 


1.492(62) 




1/2 


6 


1 -4 


4.38 


1.017(39) 


E-+ 


1/2 


3 


1 -4 


0.74 


1.365(36) 




1/2 


4 


1 - 5 


1.16 


1.357(25) 


T+- 


1/2 


8 


1 -4 


0.95 


1.276(33) 



Table 30: Results from fits to the ^ = 2, f3 = 3.3 glueball correlators on the 

10^ X 20 lattice in units of the temporal lattice spacing oj: to/ti arc used in the 
generalised eigenvalue problem, M denotes the number of operators kept after 
the truncation in C{to). 
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Channel 



to/ti M fit range x^/A^df energies 



A++ 


A 1 

91 


1 /O 

1/2 


24 


2 


7 




0.405(1^) 


A++* 


1 1 n 

i iU 


n /I 
0/1 


1 1 n 

i iU 


A 

4 


D 


U.oU 


n Ton ( 1 

U, / ZU(^10o ) 


ii ' ^ 


1 A /I 


1 /o 

1/2 


OO 

2o 


o 
2 


r 




(J. lb 


U.o75( Jo j 


1?++* 


1 A /I 

iU4 


1 /O 

1/2 


14 


1 


i) 


U.7L) 


1.183(55) 


rp+ + 


OO 


1 /9 

1/2 


1 n 

iy 


z 





U.ol 




rr++* 


06 


1 /o 

1/2 


22 


-1 
1 


o 

6 


(J. 59 


1.12o(2aJ 


A++ 
^2 


1 ^ 


1 /9 

1/2 





i 


Q 


U.UO 


l.ZlZ(4 J j 


A++* 


1 

io 


1 In 
1/2 


D 


i 


Q 

— 6 


1 1 o 

i.iz 


l.Obo^loU ) 


rp+ + 


ZO 


1 /o 
1/2 


ii 


i 


A 

— 4 


U.Uo 


l,Z4o( Jo ) 




1 ^ 

io 


1 /o 
1/2 


A 

y 


i 


q 
O 


n 1 7 
U.l / 


U, / o4( Jo j 


A h* 


15 


1/2 


9 


1 


3 


O.bl 


1.14U(112j 






1 /9 


1 /I 


i 


q 


n 97 

U.z / 






91 


1 /9 


1 1 
1 1 


1 
i 


q 


1 9n 

l.ZU 




rp- + 


( 


1/2 


o 
y 




z 


— o 


n q9 

U.oz 




rr-+* 

-^2 


/o 


1 /o 
1/2 


/? 

D 


1 

i 


Q 

o 


1 7A 


i.o75^59 J 


rp+- 


OO 


1 /o 
1/2 


1 A 
i4 


1 
i 


q 
o 


u.oy 




A A — 


11 


1/2 


5 


1 


3 


0.07 


1.224(64) 


At-* 


11 


1/2 


5 


1 


- 3 


0.79 


1.726(147) 


T+- 


22 


1/2 


9 


2 


-4 


0.33 


1.142(121) 


rr+-* 


12 


1/2 


9 


1 


- 3 


1.78 


1.663(139) 


E+- 


50 


0/1 


50 


1 


- 3 


1.01 


1.379(34) 


Tr 


80 


0/1 


80 


2 


-4 


0.57 


1.200(67) 


Tr 


80 


0/1 


80 


2 


- 5 


0.20 


1.205(81) 


A-2- 


35 


0/1 


35 


1 


- 3 


0.16 


1.352(37) 


Ar 


30 


0/1 


30 


1 


-4 


0.14 


1.703(99) 



Table 31: Results from fits to the ^ = 2, /3 = 3.5 glueball eorrelators on the 
12^ X 24 lattice in units of the temporal lattice spacing at: to/ti are used in the 
generalised eigenvalue problem, N is the number of initial operators measured 
and M denotes the number of operators kept after the truncation in C{to). 
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Cnannei 


J 


terms m the nt 


2 / Ar 




- J 1- 

A++ 





C, C2 


2.72 


4.01(15) 


A++ 





C, C2 


1.43 


7.66(71) 


E++ 


2 


C, C2 


1.06 


6.16(24) 


E++* 


2 


C, C2 


3.85 


10.59(63) 


rr++ 

^2 


2 


C, C2 


1.98 


6.14(18) 




3 


C, C2 


1.65 


11.65(57) 







c 


2.50 


6.46(18) 


E-+ 


2 


C, C2 


1.06 


8.73(34) 


rp- + 


2 


C, C2 


1.07 


8.77(39) 



Table 32: Results of the continuum extrapolations of selected glueball represen- 
tations in terms of the hadronic scale rp. The continuum spin assignment J, 
the terms in the fit, constant (c) or (ag/ro)^ (02), and the goodness of the fit, 
X^/A^DF, are also given. 





J 


terms in the fit 




rriYPc /mrp++ 




E++ 


2 


c 


1.20 


1.01(2) 


6.20(30) 


E++* 


2 


C, C2 


2.20 


1.74(12) 


10.68(105) 




2 


C, C2 


3.01 


1.70(8) 


10.44(80) 




3 


C, C2 


1.11 


1.90(11) 


11.66(102) 


T++ 


3 


C, C2 


2.60 


1.90(8) 


11.66(83) 







c 


1.39 


1.09(4) 


6.69(44) 




2 


C, C2 


1.23 


1.45(8) 


8.90(75) 


T+- 


1 


C, C2 


2.29 


1.54(7) 


9.45(71) 



Table 33: Results of the continuum extrapolations of selected glueball ratios, 
mG/mj,++. The continuum spin assignment J, the terms in the fit, constant (c) 

or {us/ro)^ (C2), and the goodness of the fit, x^/A^df, are also given. The last 
column lists the masses converted to units of r^^ using the continuum result for 
the mass of the tensor glueball T^^, values in bold face will be used further on. 
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M 


fit range 






3.0 1 


5 


3 10 


0.317(5) 


1.34 


2 


3 


4-10 


0.375(8) 


0.31 


4 


3 


3-11 


0.487(7) 


0.65 


5 


5 


4-11 


0.520(11) 


0.35 


8 


5 


3-7 


0.663(17) 


0.09 


9 


3 


4-11 


0.754(37) 


0.87 


10 


5 


4-11 


0.678(32) 


1.56 



Table 34: Collection of results of the torelon measurement at /? = 3.0, using 
the C = 4 perfect action. For each momentum = + Py we list the number 
of operators M kept after the first truncation in the variational method, the 
plateau region on which the fit of the correlators to the form Z{p'^) cxp{—tE(p'^)) 
is performed (fit range), as well as the extracted energy together with the 
per degree of freedom (x^/A^df)- 
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